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1 .0  INTRODUCTION 


This  report  summarizes  the  methods  selected  for  use  in  Missile  DATCOM. 
Since  the  preliminary  design  of  missiles  requires  the  use  of  rapid  and  accu¬ 
rate  methods,  it  became  clear  from  the  Feasibility  Study  (Reference  1)  that 
extremely  sophisticated  paneling  codes  would  be  incompatible  with  the  design 
environment.  Although  codes  exist  which  provide  the  quick  and  inexpensive 
tool  desired,  each  has  major  limitations  which  reduce  design  productivity. 

Every  effort  was  made  to  produce  a  method  collection  which  encompasses  a 
wide  range  of  configuration  geometries  and  flight  conditions,  however,  there 
remain  many  areas  which  cannot  be  addressed  with  a  preliminary  design  code. 
These  limitations  result  from  the  lack  of  applicable  methodology  and  are 
Identified  throughout  this  report.  The  lack  of  applicable  methods  may  be  due 
to  either  a  lack  of  any  known  accurate  methods  or  to  a  lack  of  sufficiently 
simple  methods  to  fit  the  preliminary  design  environment. 


2.0  METHOD  SELECTION  PROCEDURE 


The  methods  for  Missile  DATCOM  were  selected  in  four  stages: 

Stage  1:  Feasibility  Study  methods  screening 

Stage  2:  Development  Study  preselection 

Stage  3:  Method  incorporation  and  validation 

Stage  4:  Method  usage  in  a  Handbook  or  computer  code 

These  four  stages  permitted  a  wide  range  of  method  types  to  be  evaluated  and 
allowed  the  methods  being  evaluated  to  be  selected  based  upon  range  of  appli¬ 
cability  (Stage  1),  ease  of  use  and  automation  (Stage  2),  accuracy  when  com¬ 
pared  to  experiment  (Stage  3),  and  computational  time  (Stage  4). 

The  following  paragraphs  summarize  the  types  of  analysis  performed  during 
each  stage. 

2.1  FEASIBILITY  STUDY  SCREENING 

The  purpose  of  the  Feasibility  Study  (Reference  1)  was  to: 

(a)  Determine  the  range  of  geometric  and  flight  conditions  for  which 
methodology  should  be  specified, 

(b)  Survey  the  literature  for  applicable  methodology  and  determine  those 
areas  lacking  in  appropriate  methods  and  needing  further  development, 

(c)  Det.rmine  the  structure  of  a  Handbook  and/or  computer  program  which 
allows  rapid  and  accurate  use  of  the  available  methodology  in  the 
missile  design  environment,  and 

(d)  Assess  the  feasibility,  or  probability  of  success,  of  developing  a 
Missile  DATCOM  and  recommend  a  Handbook  and/or  computer  program 
format. 
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The  ranges  of  geometric  and  flight  conditions  were  determined  from  analy¬ 
sis  of  existing  experimental  data  and  existing  operational  missiles  world¬ 
wide.  The  results,  shown  in  Figure  1,  were  categorized  by  Priority.  Priority 
1  were  the  ranges  of  which  at  least  75  percent  of  the  configurations  and  their 
corresponding  flight  conditions  would  be  placed;  Priority  2  encompasses  more 
than  90  percent  of  the  total  configurations  studied. 

The  available  methods  were  assessed  using  the  eight  criteria  given  in 
Figure  2.  Generally,  methods  were  sought  which  wero: 

(a)  Theoretical  or  semi -empirical ,  so  that  the  more  "unusual"  geometries 
could  be  accurately  evaluated  and,  hence,  limit  empirical  data  base 
extrapolations, 

(b)  Require  minimum  geometry  detail, 

(c)  Demonstrated  to  be  accurate  for  the  design  environment, 

(d)  Currently  used  in  Industry  and  thought  to  be  good,  reliable  tech¬ 
niques, 

(e)  Applicable  to  a  wide  range  of  geometries  and  flight  conditions,  and 

(f)  Generally  traceable  in  origin  and  compatible  with  other  selected 
techniques. 

After  studying  more  than  315  individual  methods,  it  became  clear  that  com¬ 
ponent  build-up  methods  were  better  suited  to  fulfilling  the  eight  selection 
criteria.  It  was  also  demonstrated  that  in  spite  of  deficiencies  in  some 
areas,  sufficient  "good"  methodology  existed  to  develop  a  design  tool  appli¬ 
cable  to  preliminary  design. 

A  qualitative  assessment  of  the  two  approaches  (component  build-up  and 
panel  methods)  is  shown  in  Figure  3.  In  general,  the  cost  of  using  component 


4 


PARAMETER 

SYMBOL 

PRIORITY  1 

PRIORITY  2 

SOURCE 

ANGLE  OF  ATTACK,  DEG. 

at 

-204  <*430 

-  180<<J<180 

AEROMECHANICS 

ANGLE  OF  YAW,  DEG. 

-2040620 

-18O<0<T8O 

SURVEY  AND 

AERODYNAMIC  ROLL,  DEG. 

04|<£|445 

04|<£|4l80 

WORLD'S 

MACH  NUMBER 

M 

04M46 

04M410 

MISSILE 

BODY  FINENESS  RATIO 

d/d}8 

64(l/d)B^20 

M  !/d)j430 

SYSTEMS 

NOSE  FINENESS  RATIO 

(,/d'N 

.S^(  l/d;  -5 

C^(l/d)Ni7 

FIN  EXPOSED  SPAN  TO  DIAMETER 

b/d 

L*b/d46 

0£b/d4l0 

FIN  ASPECT  RATIO 

AR 

0.64AR44 

0.1<AR4l0 

FIN  PLANFORM 

TRIANGULAR 

ALL 

TRAPEZOIDAL  . 

WING/TAIL  ORIENTATION 

IN-LINE 

ALL 

CONTROL  METHOD 

ALL  MOVEABLE 

ALL 

FIN 

REYNOLDS  NUMBER/FT 

rn 

3xI05^Rn42x107 

103^R^3xl07 

MACH-ALTITUDE 

FIN  DEFLECTION/INCIDENCE,  DEG. 

& 

0  <8430 

0^60 

BOUNDARY 

ROLL  RATE,  RAD/SEC. 

P 

0<|  p  |=4l 

.  0^|  p[i48 

MISSILE  SYSTEM 

PITCH  RATE,  RAD/ SEC. 

q 

Oi|q|^1.5 

O^i  q  |^3 

ANALYSIS 

YAW  RATE,  RAD/SEC. 

r 

0^|  r  |^1. 5 

0^|  r  \£3 

FIN  DEFLECTION  RATE,  RAD/SEC. 

8 

o^iai^io 
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FIGURE  1  RANGE  OF  GEOMETRIC/FLIGHT  CONDITIONS 
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V  MCTHOD  APPROACH 
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-  Existing  Validation 

-  Data  Base  Available 


FIGURE  2  ASSESSMENT  CRITERIA 


FIGURE  3  APPROACHES  TO  AERODYNAMIC  PREDICTION 


build-up  methods  are  far  less  than  panel  type  methods,  due  to  the  increase  in 
required  configuration  detail  for  panel  methods.  Although  computer  costs  can 
be  a  factor  in  selection,  the  cost  in  labor  hours  is  far  more  significant; 
paneling  codes  can  require  as  much  as  ten  times  more  labor  for  the  same  con¬ 
figuration  shape. 

Paneling-type  methods  are  better  suited  for  more  arbitrary  configura¬ 
tions.  The  majority  of  component  build-up  methods  have  been  developed  for 
conventional  missile  shapes  (circular  or  elliptical  bodies  with  straight- 
tapered  fins,  and  planar  or  cruciform  fin  arrangements.  Component  build-up 
methods  are  fairly  accurate  for  conventional  configurations  and  easily  lend 
themselves  to  parametric  analysis.  For  these  reasons  and  because  most  near 
term  missiles  will  continue  to  have  conventional  shapes,  component  build-up 
methods  are  the  best  choice.  Another  plus  Is  that  with  the  build-up  approach 
It  is  easier  to  substitute  experimental  data  which  allows  for  more  detailed 
analysis. 

The  development  of  Missile  DATCOM  was  found  to  be  feasible  and  methodology 
available  to  perform  preliminary  or  conceptual  misrile  design.  Ho  comprehen¬ 
sive  collection  of  missile  design  methods  is  available,  and  the  assembly  of  a 
methods  Handbook  would  fill  this  important  void. 


Existing  computer  programs  were  available,  but  were  often  limited  In  ap¬ 
plicability,  poorly  programmed ,  too  complex,  or  substantially  undocumented. 
To  overcome  these  problems,  Missile  DATCOM  was  to  be  developed  as  follows: 


(a)  A  Handbook  and  its  companion  computer  program  were  to  be  developed  in 


i  1  *  s 

Moro i  ici 


t u  „  .....u  i , . , j . 

me  naiiuuuu*.  wuuiu  muuuc. 


/  ■»  \ 
v  * ; 


—  U  ki  i-f  4-  k  itomink  r!  {  r  #*  l  ■  r 

a  ui  ici  uuw  uiui  wuyu  u  i  Jkua 


sion  of  the  physical  phenomena  being  modeled,  (2)  a  description  of  the 
selected  method,  including  all  equations,  tables,  arid  charts  necessary  for 
its  use,  and  (3)  a  bibliography. 
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(b)  The  computer  program,  developed  concurrently  with  the  Handbook,  would 
reflect,  as  a  minimum,  the  same  capability  as  the  Handbook.  The  program 
was  to  be  written  with  the  following  characteristics:  (1)  to  conform  with 
American  National  Standards  Institute  (ANSI)  standard  FORTRAN  language  for 
ease  of  use  on  a  wide  variety  of  computer  systems,  (2)  structured  code 
using  the  concept  of  structured  programming  for  code  readability,  (3) 
liberal  Internal  documentation  to  aid  code  development  and  increase  its 
utility,  (4)  modular  method  subroutines,  so  that  they  can  be  easily  ex¬ 
changed  for  better  techniques  in  the  future  without  the  need  for  extensive 
code  revisions,  (5)  Inputs  were  to  be  structured  to  minimize  the  number 
and  type  of  inputs  required,  (6)  designed  to  minimize  computer  execution 
time  and  memory,  and  (7)  a  User's  Manual  which  defines  the  inputs  required 
and  the  outputs  available. 

2.2  DEVELOPMENT  STUDY  PRESELECTION 

In  many  cases,  the  Feasibility  Study  showed  that  more  than  one  technique 
was  available  to  compute  a  particular  aerodynamic  parameter.  Hence,  it  was 
necessary  to  select  one  for  use.  One  way  of  accomplishing  this  task  was  the 
use  of  a  methods  development  tree  where  the  assumptions  made  In  deriving  the 
techniques  could  be  directly  compared,  and  a  method  with  the  least  simplify¬ 
ing  assumptions  could  be  given  more  weight  In  the  selection  process. 

Preselection  also  considered  the  currently  available  status  of  the 
methods.  Since  the  primary  analysis  tool  was  expected  to  be  the  computer 
code,  the  availability  of  a  well  coded  and  documented  procedure  was  important 
for  overall  development  efficiency  and  future  code  maintenance. 

2.3  METHOD  INCORPORATION  AND  VALIDATION 

This  third  stage  of  method  selection  is  meant  to  automate  the  selected 
methods,  verify  they  provide  accurate  results,  and  then  compare  their  accuracy 
to  experimental  data  or  theory,  whichever  is  most  appropriate.  In  many  cases 


9 


j  the  selected  methods  already  existed  In  automated  form,  making  conversion  into 

I  the  Missile  OATCOM  code  simply  one  of  interfacing  the  inputs  and  outputs.  In 

other  cases,  new  code  had  to  be  written. 

I  Whenever  computer  code  is  used  to  represent  an  analytical  technique,  it 

i 

j  typically  takes  months  of  testing  to  reduce  programming  error.  In  order  to 

j;  shorten  the  time  needed  to  make  this  code  operational,  a  "Top-Down  Design" 

and  "Structured  Programming"  technique  was  used.  In  addition,  to  make  the 
code  as  transportable  between  machines  as  possible,  ANSI  standard  FORTRAN-66 
'  language  (ABSI  X38, 01966)  was  selected. 

j ' 

j  The  concept  of  "Top-Down  Design"  means  that  the  code  is  designed  in  a  tree 

structure,  where  the  main  program  is  the  root,  and  each  task  is  assigned  to  a 
separate  branch.  Each  branch,  or  task,  is  independent  of  others;  that  is,  the 
type  of  method  incorporated  has  no  bearing  on  the  type  of  method  selected  for 
other  tasks  in  the  system.  Hence,  individual  methods  can  be  replaced  at  will 
■  without  influencing  the  remainder  of  the  code.  The  basic  tree  structure  of 

the  code  is  illustrated  in  Figure  4.  A  significant  by-product  of  this  tree 
structure  Is  well  defined  tasks  which  can  be  overlayed  to  reduce  the  needed 
computer  core  required  for  execution. 

j 

i  I 

I I  The  second  concept,  Structured  Programming,  defines  the  way  the  code  Is  to 

>  j 

be  written,  with  emphasis  on  internal  documentation.  A  large  number  of 
i!  romment  cards  are  incorporated  in  the  program  to  determine  the  logic  and  func 

|  ional  processes. 

1 

i 

;  Both  concepts  resulted  in  the  successful  development  of  the  code.  Over 

j’|  40,000  lines  of  code  were  Implemented  in  less  than  three  years,  with  few  inain- 

!  tenance  tasks  identified  As  mentioned  earlier,  a  method  must  be  verified  as 

being  coded  correctly  before  Its  accu  ay  can  be  correctly  determined.  This 
was  performed  by  parametrically  varying  the  method  Subroutine  inputs  and  plot¬ 
ting  the  outputs  as  a  function  of  the  parametric  variable.  By  comparing  these 
results  with  theory,  or  experiment,  each  method  was  verified  as  being  an  ana¬ 
lytical  model  of  the  method  selected. 
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FIGURE  4  TOP  LEVEL  EXECUTIVE  OVERLAY  STRUCTURE 


To  evaluate  the  accuracy  of  the  individual  methods  an  accepted  accuracy 
criteria  was  employed.  The  accuracy  criteria  selected  was  the  one  developed 
by  Krieger  and  Williams  (Reference  2).  It  is  based  on  the  idea  that  the  ac¬ 
curacy  to  which  an  aerodynamic  coefficient  needs  to  be  predicted  should  be  a 
function  of  what  the  coefficient  is  going  to  be  used  to  compute.  If  it  is 
used  to  compute  obtainable  load  factor  then  it  may  not  be  required  to  be  as 
accurate  as  If  it  was  being  used  to  compute  the  range.  Therefore,  the  uses  to 
be  made  of  the  different  aerodynamic  coefficients  must  be  determined  follow¬ 
ed  by  a  determination  of  accuracy  required  for  the  performance  parameters. 
The  accuracy  required  In  the  prediction  of  the  aerodynamic  coefficients  is 
then  backed  out  of  the  equations  for  the  performance  parameters.  This  process 
is  outlined  for  the  longitudinal  coefficients  C^,  and  in  Figure 

5.  From  this  accuracy  criteria  it  is  found  that  it  is  acceptable  to  predict 
C,,  and  C  within  +20  percent  while  C.  is  often  required  to  be  predicted 
to  within  +10  percent.  Application  of  this  method  allowed  for  a  quantitative 
assessment  of  the  accuracy  of  the  Individual  methods. 

In  addition  to  individual  method  accuracy,  it  was  necessary  to  insure  that 
they  were  compatible  as  a  function  of  Mach  number.  As  illustrated  in  Figure 

6,  It  was  desirable  that  there  be  no  predictive  discontinuities  at  the  limits 
of  method  applicability.  Even  though  this  was  not  practical,  the  differences 
at  method  switchover  were  kept  at  a  minimum.  Method  switchover  occurred  at 
the  following  conditions: 


CONFIGURATION  MACH  NUMBER  COEFFICIENTS 

BODY  1.2  ALL 

2.0  ALL 

FINS  0.8  ALL 

1.4  ALL 


The  reasons  for  the  switchover  are  defined  in  the  method  selections  described 
in  Section  3  and  Section  4. 
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FIGURE  5  DERIVATION  OF  ACCURACY  CRITERIA 


CHECKS  METHODS  COMPATIBILITY 


2.4  COD£  AND  HANDBOOK  METHOD  SELECTION 


The  goals  for  the  Missile  DATCOM  program  were  to  have  the  computer  code 
and  the  Handbook  complement  each  other,  yet  be  independent  so  that  either 
could  be  used  without  requiring  the  use  of  the  other.  In  some  cases,  the 
"best"  methodology  could  be  efficiently  used  on  computer  equipment  and  not 
be  easily  solved  by  hand  in  a  methods  Handbook.  It  w.»uld  have  been  possible 
to  use  the  code  to  generate  Handbook  design  charts,  which  was  ultimately  done 
for  the  skin  friction  calculation  method,  but  this  approach  was  avoided  since 
It  does  not  lend  Itself  to  a  thorough  understanding  of  the  method.  By  avoid¬ 
ing  this,  the  Handbook  and  code  methods  are  Independent  checks  of  one  another, 
and  In  some  instances  are  from  different  method  sources,  so  the  methods  in  the 
Handbook  could  be  kept  as  simple  as  possible. 

Where  there  are  method  differences  between  the  Handbook  and  the  computer 
code,  those  differences  are  described  in  the  method  description  sections  of 
this  report. 
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3.0  BODY  ALONE  METHOD  SELECTIONS 


This  section  summarizes  the  methods  selected  for  body  alone  aerodynamics. 
It  is  divided  into  two  sections:  bodies  with  circular  cross  sections  (axisym- 
metric  bodies)  and  bodies  with  noncircular  cross  sections.  Each  section  is 
further  subdivided  by  aerodynamic  coefficient. 

It  should  be  noted  that  a  significant  number  of  "easy  to  apply"  methods 
for  noncircular  bodies  have  been  used  for  bodies  with  elliptically-shaped 
cross  sections.  Although  methods  can  be  extended  to  handle  more  arbitrary 
shapes,  this  has  not  been  done  due  to  the  limited  amount  of  information  avail¬ 
able  for  extending  these  techniques.  The  method  used  in  the  computer  code 
for  bodies  with  arbitrary  cross  sections  is  numerically  intensive,  therefore 
no  hand  calculation  methods  are  presented. 

3.1  BODIES  WITH  CIRCULAR  CROSS  SECTIONS 

This  section  summarizes  the  body  alone  methods  investigated  for  Inclusion 
into  the  Missile  OATCOM  methods  compendium.  Approximately  108  individual 
methods  were  Identified  during  the  Missile  DATCOM  Feasibility  Study  as  either 
recommendations  or  alternates 

Method  selection  began  with  the  Feasibility  Study  results.  Each  of  the 
recommended  techniques  were  again  studied  using  the  assessment  criteria  shown 
In  Figure  2.  The  final  selections  are  given  in  Figure  7,  and  a  summary  com¬ 
parison  with  the  Feasibility  Study  recommendations  Is  given  in  Table  1.  In 
nearly  all  cases,  the  Feasibility  Study  recommendations  were  the  final  method 
selections. 

As  shown  in  Figure  8,  component  build-up  methods  remain  questionable  for 
some  configurations.  For  example,  the  subsonic/transonic  axisymmetric  methods 
included  are  only  applicable  to  conical  or  tangent  ogive  nose  shapes  attached 
to  cylindrical  bodies,  whose  lengths  are  at  least  six  body  diameters.  The 
reasons  for  the  configuration  limitations  are  as  follows: 
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FEASIBILITY  RECOMMENDATION  METHOD  SELECTED  REASON  FOR  CHANGE 
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1  TANGENT  OGIVES:  y/ 

FIGURE  8  COMPONENT  BUILD-UP  APPLICABILITY 

(a)  In  the  subsonic  and  transonic  Mach  regimes  no  easy  to  apply  method  Is 
available  to  analyze  general  nose  shape  contours.  Although  design  charts 
exist  for  selected  nose  shapes,  the  methods  require  substantial  input  de¬ 
scription  and  user  Intervention.  Therefore  they  are  not  Included  In  the 
Missile  DATCOM  code.  Parametric,  results  in  the  form  of  design  charts  have 
not  been  published  and  are  therefore  unavailable  for  Incorporation. 

(h)  Methods  to  compute  the  effect  of  spherical  nose  bluntness  Is  only 
available  for  CA  at  zero  angle  of  attack.  The  blunted  cone  method  used 
Is  weak.  (Design  charts  are  available  for  blunted  cones,  but  are  not  In 
the  public  domain.) 


(c)  Elllptlcally  shaped  nose  bluntness  was  not  addressed  in  any  litera 
ture  sources  searched. 


(d)  The  extreme  dependency  of  the  nose  flow  field  on  the  afterbody  for 
short  cylindrical  afterbodies  prohibits  the  generalization  of  a  method  for 
arbitrary  afterbody  shape  or  nose-flare  geometries. 

(e)  The  existing  flare  methodology  at  subsonic  and  transonic  speeds  is 
limited  to  conical  flares  whose  surface  slope  angle  Is  small. 

The  following  paragraphs  detail  each  of  the  aerodynamic  coefficient 
methods  selected  and  their  limitations. 

3.1.1  Axial  Force  Coefficient 

In  general,  the  methods  recommended  during  the  Feasibility  Study  were 
found  to  be  acceptable  during  method  selection.  Each  of  the  following  para¬ 
graphs  summarize  the  results  observed  during  the  method  Implementation  and 
validation. 

3. 1.1.1  Skin  Friction  -  The  method  of  Van  Driest  (referred  to  as  Van  Driest 
IT)  was  found  to  produce  acceptable  skin  friction  drag  characteristics  for 
turbulent  boundary  layers  at  all  speeds  (Reference  3),  The  method  of  Blaslus, 
for  laminar  boundary  layers,  and  the  MDAC  method  for  laminar  boundary  layer 
transition  completed  the  scope  of  the  skin  friction  calculations  (References  4 
and  5  respectively).  Since  skin  friction  is  difficult  to  extract  from  a  ma¬ 
jority  of  experimental  results  for  bodies,  the  absolute  magnitude  of  the 
method's  predictive  accuracy  cannot  be  determined.  However,  based  upon  total 
drag  comparisons  with  experimental  data,  it  was  deduced  that  the  skin  fric¬ 
tion  calculations  are  quite  adeqi  ite  for  preliminary  design  purposes. 

Since  the  Van  Driest  II  method  is  iterative,  it  cannot  be  easily  computed 
by  hand  using  a  methods  handbook.  However,  since  it  includes  the  effects  of 
wall  temperature,  Mach  number,  Reynolds  number,  and  surface  roughness  It  far 
surpasses  the  available  simple  techniques.  Hence,  the  production  of  design 
charts  for  the  methods  handbook  was  considered  the  best  compromise;  an  example 
Is  shown  in  Figure  9.  In  the  computer  code  the  free  stream  temperature  is 
used  as  the  wall  temperature. 
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REYNOLDS  NUMBER 

FIGURE  9  SKIN  FRICTION  METHOD  RESULTS 


3 . 1 . 1 . 2  Subsonic  Pressure  Drill  -  The  Method  of  Hoerner,  which  correlates 
pressure  drag  to  skin  friction  drag  was  retained  from  the  Feasibility  Study 
results  (Reference  6).  There  were  no  subsonic  (Mach  numbers  less  than  0.6) 
predictive  comparisons  which  demonstrated  that  this  approach  is  inappropriate. 


3.1 ,1.3  Transonic  Pressure/Wave  Drag  -  For  ogive-cylinder  configurations  the 
method  of  Chaussee  (Reference  7)  was  initially  selected  because  of  its  basis 
In  transonic  flow  theory.  However,  after  the  Feasibility  Study  was  completed, 
It  was  determined  that  the  Chaussee  results  were  substantially  In  error  for 
selected  nose  fineness  ratio  combinations.  These  discrepancies  were  corrected 
by  Devan  (Reference  8)  and  form  the  basis  of  the  final  method  selected.  Com¬ 
parison  with  experiment  has  shown  that  this  technique  gives  excellent  results 
(Figure  10),  although  It  was  required  to  be  coded  as  a  multivariable  table 
lookup.  The  original  Lagrange  Interpolation  was  replaced  with  a  linear  inter¬ 


polation  to  reduce  Interpolation  error. 
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FIGURE  10  SPHERICALLY  BLUfITED  HOSE  DRAG  RESULTS 
USING  DE VAN'S  METHOD 


The  recommended  approach  for  cone-cylinder  pressure/wave  dreg  was  develop¬ 
ed  by  MDAC.  It  covered  a  wide  range  of  cone  fineness  ratios  and  nose  blunt¬ 
ness  ratios.  Although  this  technique  Is  probably  one  of  the  better  approach¬ 
es,  the  method  Is  considered  propletary  and  was  not  used.  The  alternative 
approach  was  one  derived  by  Devan  from  experimental  observations,  which  Indi¬ 
cated  a  very  small  effect  of  nose  bluntness.  Since  this  approach  is  not  en¬ 
tirely  accurate,  nor  does  it  cover  a  sufficiently  wide  range  of  nose  fineness 
ratios,  It  is  recommended  that  it  be  replaced  with  a  better  technique  when 
available.  The  coded  method's  limitations  are  summarized  In  Figure  11. 


3. 1.1. 4  Supersonic  Wave  Drag  -  Two  techniques  were  selected  for  body  wave 
drag,  the  nullified  Second  Order  Shock  Expansion  (SOSE)  method  of  DeuarPiette 
(Reference  9),  and  the  Van  Dyke  hybrid  theory  method  as  automated  by  NSWC  for 
the  Aeropredlctlon  Code  (Reference  10).  Both  techniques  were  selected  so  that 
the  entire  supersonic  Mach  regime  could  be  adequately  covered. 
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FIGURE  11  NOSE-CYLINDER  GEOMETRY  LIMITATIONS 
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At  Mach  numbers  below  approximately  2,0,  there  is  a  tendency  for  SOSE  to 
overpredict  wav  drag  by  approximately  20  percent  (See  Figure  12).  Since  this 
approach  models  geometry  as  a  series  of  cone  frastrums,  this  type  of  error  can 
be  expected  at  low  supersonic  speeds  for  some  shapes.  At  a  Mach  number  of  ap¬ 
proximately  1.4  or  less,  the  SOSE  code  will  compute  subsonic  flow  behind  the 
shock,  with  large  nose  tip  pressures  being  a  result;  hence,  the  axial  force 
will  climb  to  extremely  high  values.  The  flow  expands  to  supersonic  speeds 

with  no  noticeable  effect  on  C..  and  C  curve  slope  near  zero  angle  of 

N  JTI 

attack. 


To  avoid  this  problem  at  low  supersonic  speeds,  the  Van  Dyke  Hybrid  method 
was  Incorporated.  Its  use  eliminated  the  large  wave  drag  characteristic  of 
SOSE,  but  at  the  expense  of  computational  time.  Although  Hybrid  Theory  im¬ 
proved  the  prediction  accuracy  at  lower  speeds,  it  does  not  compute  the 
effects  of  spherical  nose  bluntness  accurately,  and  for  this  reason  is  not 
recommended  for  blunted  noses  at  any  speed.  The  effect  of  nose  bluntness  at 
low  supersonic  speeds  is  a  noticable  methods  problem  area.  It  is  recommended 
that  a  new  method  be  Incorporated  in  the  Mach  range  from  1.2  to  1.5  for 
blunted  noses.  The  easiest  approach  would  be  to  extend  the  transonic  methods 
from  the  current  Mach  1.2  to  Mach  1.5;  sufficient  data  exists  In  the  transonic 
method  data  tables  to  permit  this  extension. 


3. 1.1. 5  Boattail  Pressure/Wave  Drag  -  The  Payne  correlation  of  boattail  drag 
at  subsonic  and  transonic  speeds  is  comprehensive  and  is  restricted  only  by 
the  limitations  detailed  in  Figure  13  (Reference  11).  It  was  found  to  be  suf¬ 
ficiently  accurate  for  preliminary  design.  However,  It  does  assume  that  suf¬ 
ficient  distance  exists  between  the  nose  and  the  boattail  to  minimize  nose 
flow  field  coupling  (approximately  six  calibers).  Shorter  cylindrical  after¬ 
bodies  should  not  be  investigated  using  this  technique,  and  no  other  method  Is 
available  that  can  be  substituted.  Hence,  the  derivation  of  a  new  subsonic/ 
transonic  boattail  drag  model  is  required.  A  recommended  approach  would  be  to 
use  the  code  RAXBQB  to  generate  parametric  results  which  could  be  interpolated 
for  the  computer  code,  and  used  to  create  a  design  chart  for  the  methods  hand¬ 
book  (Reference  12). 
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FIGURE  13  BOATTAIL  GEOMETRY  LIMFATIOKS 

Comparisons  with  experimental  data  at  supersonic  speeds  have  shown  that 
the  SOSE  and/or  Hybrid  methods  are  sufficiently  accurate  for  preliminary  de¬ 
sign  since  they  Inherently  Include  nose/afterbody  coupling. 

3. 1.1. 6  Flare  Pressure/Wave  Drag  -  The  Feasibility  Study  did  not  Ident¬ 
ify  a  method  for  flares  In  the  subsonic/transonic  regime  slnca  It  was  thought 
that  their  usage  was  limited  to  projectile  shapes.  Subsequent  Information 
revealed  that  there  exists  much  more  Interest,  than  first  believed,  so  two 
methods  were  later  selected  for  stuty.  The  first  was  an  empirical  data  base 
of  a  wide  variety  of  flare  fineness  ratios;  It  was  rejected  because  of  the 
extremely  short  cylindrical  centerbody  separating  It  from  the  nose.  The 
second,  a  U.S.  Army  correlation,  seemed  to  be  a  reasonable  alternative  al¬ 
though  It  only  addressed  conical  flares  (Reference  13).  The  limitations 
placed  on  the  geometry  are  described  in  Figure  14.  Since  no  Other  method 
exists  it  was  selected.  It  Is  recommended  that  an  improved  new  technique  be 
developed. 
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FIGURE  14  FLARE  GEOMETRY  LIMITATIONS 


At  supersonic  speeds,  flares  can  be  analyzed  quite  easily  using  either  the 
SOSE  or  Hybrid  methods  described  above.  Comparison  with  experiment  has  shown 
the  pressure  distributions  calculated  to  be  quite  accurate  at  lower  angles  of 
attack  (less  than  eight,  degrees). 


3. 1.1. 7  General  Nose  Shape  Pressure/Have  Drag  -  The  only  nose  shapes  cur¬ 
rently  being  addressed  at  subsonic/transonic  speeds  are  tangent  ogives  and 
cones.  There  Is  ar  urgent  need  to  develop  methods  for  other  shapes,  such  as 
Ksack  series,  power-law,  and  secant  ogive  noses.  The  approach  recommended 
Is  to  use  the  code  RAXBO  and  parametrically  vary  the  nose  fineness  ratio  to 
develop  design  charts.  The  supersonic  speed  regime  Is  once  again  adequately 
covered  by  the  SOSE/Hybrld  combination. 


29 


3. 1.1. 8  Base  Drag  -  The  method  selected  by  the  Feasibility  Report  was 
slightly  changed  simply  due  to  automation  convenience.  The  use  of  NASA-TR-R- 
100  Is  thought  to  be  slightly  more  advantageous  due  to  the  wide  range  of 
flight  test  results  used  In  derivation  of  the  method  (Reference  14).  The 
method  was  modified  using  a  U.  S.  Army  method  so  that  the  base  drag  due  to 
boattalls  and  flares  could  be  accommodated  (Reference  13). 

3. 1.1. 9  Jet  Exhaust  and  Protuberance  Effects  -  Although  methods  were  Ident¬ 
ified  In  the  Feasibility  Study  for  the  effects  of  jet  exhaust  and  protuber¬ 
ances,  they  proved  to  be  too  limited  In  application  or  too  general  to  be 
useful.  Hence,  neither  set  of  techniques  were  Incorporated. 

3.1.1.10  Anole-of-Attack  Effect  -  The  recommended  Allen  and  Perkins  plus  Jor¬ 
gensen  method  was  Incorporated  (Reference  IS).  Prediction  comparisons  showed 
that  the  techniques  were  adequate  to  approximately  30  degree  angle  of  attack. 

However,  In  the  transonic  Mach  regime  the  discrepancy  with  experiment  was 
largest  between  15  degrees  and  60  degrees  angle  of  attack.  Some  predictions 
followed  the  experimental  results  quite  accurately,  while  others  were  signi¬ 
ficantly  In  error.  It  is  not  known  whether  this  effect  Is  Reynolds  number, 
nose  shape,  or  total  body  fineness  ratio  dependent.  Analyses  which  are  more 
dependent  upon  results  In  this  angle  of  Attack  range  will  require  a  signifi¬ 
cantly  improved  axial  force  due  to  angle  of  attack  model. 


3.1.2  Normal  Force  and  Pitching  Woment  Coefficient 


The  methods  selected  for  configuration  normal  force  and  pitching  moment 

_ _ r _  _ A  ~  _  X.  XL  _ _ _ - - J  i.  XL.  r _ JLi  lJx..  . 

are  lar  bUpvnur  iu  uivse  uriymaiiy  r  otuininiiuuu  m  trie  rvaaiui'ity  otuuy. 

This  is  due  to  the  discovery  of  methods  not  previously  available  In  the  open 
literature  and  the  Improvement  of  the  selected  supersonic  techniques. 
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It  Is  assumed  that  the  normal  force  and  pitching  moment  coefficients  can 
be  approximated  as  the  sum  of  potential  and  viscous  contributions: 

CN  '  CNP  *  CNV 
"  ^rap  +  CmV 

where, 

Cu  -  0.5  Cu  sln(So)  cos  (a/2) 

Np  Na 

CNV  "  n  Cdc  SP  sin  ^/Sref 

C  -  0.5  C  s1n(2a)  cos (a/2) 
mp  ma 

C»v  “  n  Cdc  SP  sin  ^a^Sref*Xcg~XC*Lref 

3. 1.2.1  Potential  normal  Force  and  Pitching  Moment  -  Two  separate  methods 
wer  Implemented:  the  Second  Order  Shock  Expansion  (50SE)  method  and  the  Van 
Dyke  Hybrid  Theory.  As  Illustrated  in  Figure  15,  a  significant  methods 
overlap  region  exists  for  the  most  commonly  used  noses.  It  would  appear  ad¬ 
vantageous  to  Include  only  the  SOSE  method  at  all  supersonic  speeds  and  sim¬ 
plify  the  code.  At  Mach  2.3,  the  SOSE  pressure  coefficient  distributions  over 
a  3:1  tangent  ogive  at  eight  dcgreos  eng le-of -attack,  compare  favorably  to  the 
more  complex  Hybrid  and  Shock  Capturing  codes.  Figure  16.  Other  SOSE  compari¬ 
sons  to  the  configurations  shown  In  Figure  17  were  also  quite  satisfactory 
(Figures  18  through  24). 

Using  a  large  data  base  of  experimental  results,  the  accuracy  of  the 
SOSE  method  was  assessed.  The  results,  shown  in  Figures  25  and  26,  divided  by 
nose  shape,  shows  that  a  large  percentage  of  the  configurations  (over  65  per¬ 
cent)  can  be  predicted  with  10  percent  error  or  less,  and  over  75  percent  can 
be  predicted  with  less  than  15X  error.  The  average  error  on  normal  force  and 
pitching  moment  coefficient  slope,  summarized  In  Figure  27,  is  far  below  the  ± 
20  percent  allowable  using  the  accuracy  criteria  developed  In  Section  2. 
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FIGURE  15  RECURS  OF  APPUCABIUTY-^ifERSONIC  YMRQRIfS 

In  Figure  28  these  results  were  divided  by  the  Most  cossswn  geoMBtry  cri¬ 
teria 


(1)  e/fN<0.4 

(2)  H  <  2 

(3)  f A  -  0 

(4)  M  >  2,  fA  >  0 


-  an  assumed  poor  SOSE  application  region, 

-  an  assumed  lower  Mach  limit  for  SOSE, 

-  no  cylinder  added  to  the  nose,  and 

-  the  best  SOSE  application  region. 


As  assumed,  criteria  4  was  the  best  demonstrated  operating  region  for 
SOSE.  However,  the  average  error  for  the  other  regions  was  a  scant  2.0  per¬ 
cent  to  3.4  percent  greater. 


FIGURE  16  MISSILE  DATCOM  COMPUTED  PRESSURE  COMPARISONS 
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LONGITUDINAL  STATION 
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MACH- 1.6  ALPHA-8. 0  PHI 


ALPHA-8.0  PHI-0. 


LONGITUDINAL  STATION 

GURE  ?2  SOSE  PRESSURE  COEFFICIENTS  FOR  A  CIRCULAR-ARC-CIRCULAR-ARC 


MACH-4.63  ALPHA-8.0  PHI 


LONGITUDINAL  STATION 


ALPHA-8.0  PHI 


CONFIGURATION 

CNa  AVERAGE  ERROR 

SHARP  CONE-CYLINDER 

SHARP  TANGENT  OGIVE-CYLINDER 

9.2% 

6.2 

AVERAGE  ERROR 


6.4% 

7.6 


FIGURE  27  PREDICTION  ERROR  OF  POTENTIAL  SLOPES  USING  SOSE 


CONFIGURATION  SELECTION 

CRITERIA 

AVERAGE  ERROR 

/J/f>N<0. 4 

7.7% 

M<2.0 

7.0% 

fA=0. 

8.4% 

M>2,  fA>0. 

5.0% 

FIGURE  28  PREDICTION  ERROR  FOR  TANGENT  OGIVE-CYLINDER  CONFIGURATIONS,  CNa 
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The  same  data  base  was  used  to  demonstrate  the  accuracy  of  the  Hybrid 
Theory.  As  shown  in  Figure  29,  the  Hybrid  Theory  was  shown  to  be  approximate¬ 
ly  twice  as  accurate  as  SOSE.  Because  of  the  increased  accuracy  of  this 
method,  it  was  advantageous  to  include  it;  its  retention  also  provided  a  solu¬ 
tion  to  the  wave  drag  problem  experienced  with  SOSE  at  lower  supersonic  speeds 
(see  Figure  12) . 

The  subsonic/transonic  incremental  normal  force  slope  for  a  flare  is 
assumed  by  the  selected  method  to  be  dependent  upon  the  cylinder  diameter  to 
flare  base  diameter  ratio,  and  is  insensitive  to  a  flare  angle  for  angles  up 
to  ten  degrees  (Reference  13).  This  theory  gives  acceptable  results  for  pre¬ 
liminary  designs  as  shown  in  Figure  30.  No  systematic  trend  in  the  experi¬ 
mental  data  could  be  found  which  would  allow  a  more  accurate  method  to  be 
developed,  nor  were  data  for  higher  flare  angles  found  in  the  open  litera¬ 
ture.  Hence,  the  selected  method  appears  to  be  the  best  available  for  the 
range  of  configurations  anticipated.  However,  a  better  technique  should  be 
sought  which  includes  the  effect  of  larger  flare  angles  and  the  influence  of 
centerbody  length  since  these  are  significant  effects. 


3. 1.2.2  Viscous  Normal  Force  and  Pitching  Moment  -  The  effect  of  angle  of 
attack  is  modeled  by  assuming  that  the  drag  of  an  infinite  cylinder  normal  to 
the  flow  can  be  linearily  superimposed  upon  the  configuration  potential 
normal  force  or  pitching  moment.  Since  there  are  at  least  nine  different 
cross  flow  drag  models  available,  it  is  necessary  to  pair  the  potential  model 
with  its  corresponding  viscous  model  to  obtain  a  correct  method.  This  "pair¬ 
ing"  is  required  since  an  empirically  derived  correction  factor  is  applied  to 
the  viscous  model  in  the  subsonic/transonic  noth  regime.  The  correction 
factor,  is  the  cross  flow  drag  proportionality  factor;  It  specifies  how  two- 
dinensional  the  flow  around  the  body  becomes. 


The  choices  for  cross  flow  drag  were  reduced  to  the  three  shown  in  Figure 
31.  The  Jorgensen  model  is,  perhaps,  the  most  "empirically  correct"  of  the 
three,  since  two-dimensional  cylinder  cross  flow  data  was  used  and  not  data 
"corrected"  to  match  the  potential  model  employed.  Jorgenson  used  a 
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FIGURE  29  HYBRID  VERSUS  SOSE  PREDICTIVE  ERROR— Ctfa 
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30  SUBSONI 


METHOD  DIFFERENCES 


COMPARISON  OF  CROSS-FLOW  DRAG  MODELS  AVAILABLE 


nose-cylinder  normal  force  coefficient  slope  of  two  per  radian  to  derive  his 
cross  flow  drag  coefficient.  Moore  and  Balter  relied  on  experimental  nose  cy¬ 
linder  normal  force  coefficient  slope  to  define  the  cross  flow  drag  coeffi¬ 
cient  as  a  function  of  cross  flow  Mach  number. 

In  Figure  32,  experimental  data  for  the  cross  flow  drag  coefficient  of 
Infinte  cylinders  as  a  function  of  cross  flow  Mach  number  is  shown.  This  data 
was  then  used  to  compute  the  viscous  normal  force  and  pitching  moment  coeffi¬ 
cient.  The  dotted  line  labeled  "old  fairing*1  is  the  cross  flow  drag  coeffi¬ 
cient  variation  first  assumed.  Extensive  correlations  with  experimental  data 
showed  that  this  model  tended  to  overpredict  CN  from  10  to  20  percent  for 
typical  missile  angles  of  attack  and  Mach  number  ranges.  By  adjusting  the 
cross  flow  drag  coefficient  variation  to  the  curve  labeled  "new  fairing"  the 
error  in  most  predictions  was  held  to  10  percent  or  less.  This  modification 
Is  a  methods  calibration.  A  different  choice  of  a  normal  force  coefficient 
slope  prediction  method  would  probably  not  yield  the  accuracy  level  observed 
with  the  current  method. 


Calibrating  the  method  does  not  violate  any  method  derivation  assumptions 
or  theoretical  base.  The  assumption  that  the  normal  force  and  pitching  moment 
can  be  estimated  by  summing  potential  and  viscous  contributions  Is  based  on 
the  assumption  that  the  longitudinal  and  normal  force  vectors  can  be  vector- 
ally  summed.  This  Is  a  basic  assumption  of  linearity  for  a  physical  phenome¬ 
non  which  is  not  completely  accurate,  but  is  a  reasonable  approximation  to  a 
complex  analysis. 


H  *1 


Lateral— Di recti on si 


The  methods  described  in  Section  3.1.2  for  normal  force  and  pitching  mo¬ 
ment  are  directly  applicable  to  side  force  and  yawing  moment.  The  angle  of 
attack  used  by  those  methods  is  really  the  total  body  angle-of-att;  ck.  Any 
directional  effect  arises  from  the  axis  system  employed,  Ignoring  the  effect 
of  asymmetric  vortex  shedding.  If  CN'  and  C  1  represent  the  body  normal 
force  and  pitching  moment  coefficients,  respectively,  as  a  function  of  total 
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FIGURE  32  MODIFICATION  OF  CROSS-FLOW  DRAG 


angle-of -attack,  then  the  body  axis  normal  force,  side  force,  pitching  moment, 
and  yawing  moment  are  respectively: 


CN  -  CN'  cos(4>') 

CY  "  _CN'  sin  (♦*) 

C  -  C  '  cos(4>' ) 
m  m 

C  -  -C  '  sin(4>') 
n  m 


where. 


-1 

*  t^n  (tan  S/tan  c*) 

Methods  for  estimating  the  effect  of  "phantom  yaw"  or  asymmetric  vortex 
shedding  have  not  been  Included  Into  the  handbook  or  the  code  due  to  the  com¬ 
plexity  of  the  solution  and  the  uncertainty  of  such  flow  parameters  as  Mach 
number,  Reynolds  number,  and  angle  of  attack. 

3.2  BODIES  WITH  NONCIRCULAR  CROSS  SECTIONS 

The  methods  for  noncircular  bodies  are  variants  of  the  circular  body 
method: .  Not  only  are  these  types  of  methods  commonly  used  to  estimate  the 
aerodynamics  of  noncircular  bodies  but  they  ar8,  for  the  most  part,  the  only 
easy  to  apply  techniques  available. 

At  subsonic  speeds  one  could  model  noncircular  bodies  using  a  potential 
code,  where  the  geometry  is  modeled  as  a  matrix  of  sources  and  sinks;  3-D  Neu¬ 
mann  (Reference  16)  is  an  example  of  such  an  approach.  However,  the  costs  for 
such  a  solution  are  extremely  high,  and  the  setup  time  to  prepare  the  inputs 
are  highly  labor  Intensive.  Both  of  these  disadvantages  make  their  use  in  the 
preliminary  design  envirnoment  questionable,  and  perhaps  unnecessary  If  a 
"simple"  approach  produces  results  of  sufficient  accuracy. 
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At  transonic  speeds,  the  accurate  prediction  of  aerodynamics  for  noncir¬ 
cular  bodies  must  rely  on  very  complex  codes,  such  as  the  Implementation  of  a 
Navler-Stokes  solution.  These  types  of  analyses  are  inapproplate  for  prelim¬ 
inary  design. 

At  high  supersonic  speeds,  the  Supersonic/Hypersonic  Arbitrary  Body 
(S/HABP)  code  (Reference  17)  Is  a  reasonable  alternative  for  noncircular 
bodies  due  to  the  less  complicated  analyses  performed.  Although  this  approach 
was  selected  to  model  airbreathing  Inlets,  the  complexity  of  user  preparation 
was  extensive.  Results  using  the  Jorgensen  approach,  which  modifies  the  po¬ 
tential  and  viscous  components  of  normal  force  and  pitching  moment,  was  found 
to  be  highly  successful  for  elllptlcally  shaped  geometries,  making  the  use  of 
a  S/HABP  analysis  unnecessary.  An  extension  of  slender  body  theory  using 
apparent  area  calculations  was  found  to  be  useful  in  predicting  the  longitud¬ 
inal  aerodynamics  of  bodies  with  many  arbitrary  cross  sections.  While  the 
method  implemented  can  handle  bodies  not  normally  thought  of  as  slender  It 
does  have  serious  limitations.  These  are  discussed  In  Section  3.2.2. 

The  following  paragraphs  summarize  the  method  changes  applied  to  the  cir¬ 
cular  body  methods,  described  In  Section  3.1,  for  application  to  bodies  whose 
cross  section  can  be  represented  by  an  ellipse.  Following  the  description  of 
the  elliptic  body  methods  are  paragraphs  describing  the  arbitrary  cross  sec¬ 
tion  methods  used. 

3.2.1  Bodies  with  Elliptic  Cross  Sections 

3. 2. 1.1  Axial  Force  Coefficient  -  The  methods  selected  for  computing  the 
axial  force  of  bodies  with  noncircular  cross  sections  are  those  selected  for 
circular  bodies,  as  described  In  Section  3.1.1.  The  complexity  of  the  flow 
field  makes  accurate  solutions  extremely  costly  due  to  the  complex  analysis 
codes  required.  Hence,  the  use  of  circular  body  methods,  with  approplate  cor¬ 
rections  to  account  for  cross  section  shape,  is  the  most  economical  approach. 
The  methods  and  cross  section  correction  factors  are  described  below. 
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3. 2. 1.1.1  Skin  Friction  -  The  method  described  in  Section  3. 1.1.1  is  used 
for  noncircular  bodies.  Even  though  the  local  skin  friction  coefficient  could 
approach  flat  plate  values  over  some  portions  of  the  configuration,  the  three- 
dimensional  correction  to  flat  plate  frictlos  is  used  to  approximate  the  skin 
friction  over  the  whole  wetted  surface  of  the  geometry. 

3. 2. 1.1. 2  Subsonic  Pressure  Drag  -  Subsonic  pressure  drag  for  circular  bodies 
Is  computed  using  the  drag  due  to  skin  friction  and  a  "form  factor"  adjust¬ 
ment.  This  "form  factor"  was  empirically  derived  and  is  a  function  of  body 
fineness  ratio.  It  is  assumed  that  the  noncircular  body  fineness  ratio  (based 
upon  an  equivalent  circular  body  diameter)  can  be  used  to  estimate  the  pres¬ 
sure  drag  of  noncircular  bodies.  This  approximation  Is  adequate  for  prelim¬ 
inary  design. 

3. 2. 1.1. 3  Transonic  Pressure/Wave  Drag  -  Estimating  the  transonic  drag  char¬ 
acteristics  of  circular  bodies  Is  extremely  difficult  since  the  methods  avail¬ 
able  are  limited  to  a  small  number  of  geometries.  The  prediction  of  noncir¬ 
cular  bodies  Is  even  more  uncertain.  It  Is  assumed  that  the  pressure/wave 
drag  characteristics  would  be  the  same  as  those  for  a  circular  body  with  the 
same  longitudinal  cross-sectional  area  distribution. 

3. 2. 1.1. 4  Supersonic  Wave  Drag  -  This  is  the  only  component  of  axial  force 
for  which  a  method  specifically  Intended  for  elliptical  bodies  exists.  The 
method  derived  by  Van  Dyke  using  second-order  slender  body  theory,  does  not 
predict  a  significant  variation  In  wave  drag  for  typical  major  to  minor  axis 
ratios,  as  shown  In  Figure  33.  The  Van  Dyke  method  is  complex,  and  the  effect 
of  cross  section  shape  on  wave  drag  can  be  assumed  to  he  second  order  in  pre¬ 
liminary  design,  and  Is  Ignored  In  Missile  DATCOM's  elliptic  body  methods. 

3. 2. 1.1. 5  Boattall  Pressure/Wave  Drag  -  It  Is  assumed  that  the  Influence  of 
boattalls  are  similar  to  those  for  circular  bodies.  An  equivalent  circular 
body  has  been  assumed. 
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3. 2. 1.1. 6  Flare  Pressure/Wave  Drag  -  It  is  assumed  that  the  influence  of 
flares  are  similar  to  those  for  circular  bodies.  An  equivalent  circular  body 
has  been  assumed. 

3. 2. 1.1. 7  General  Nose  Shape  Pressure/Wave  Drag  -•  As  described  in  Section 
3. 1.1. 7,  few  easy  to  apply  methods  exist  tor  prodlcting  the  effect  of  arbi¬ 
trary  longitudinal  variation  in  the  nose  shapes  of  circular  bodies.  Being 
able  to  analyze  non-circular  bodies  is  more  difficult.  The  methods  Implement¬ 
ed  are  described  In  Section  3.2.4. 

3. 2. 1.1. 8  Base  Drag  -  Refer  to  Section  3. 1.1. 8  for  the  circular  body  methods 
since  they  are  applied  to  noncircular  geometries. 

3. 2. 1.1. 9  Angle  of  Attack  Effect  -  The  methods  for  circular  bodies,  described 
In  Section  3.1.1.10,  are  assumed  to  be  applicable  to  noncircular  geometries. 

3 .2. 1.2  Normal  Force  and  Pitching  Moment  Coefficient  -  The  circular  body 
methods  described  in  Section  3.1.2  are  assumed  to  be  applicable  to  noncir¬ 
cular  geometries,  provided  that  the  effect  of  cross  section  is  accounted  for 
using  correction  factors: 


A  slender  body  theory  (SB)  correction  is  made  to  the  potential  component 
of  each  coefficient,  and  a  Newtonian  theory  (NT)  correction  is  made  to  the 
viscous  component.  For  bodies  whose  cross  section  shape  varies  along  its 
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length,  the  constants  can  be  replaced  with  Integrals  along  the  body  length. 

It  Is  this  form  of  the  equations  which  are  referred  to  as  the  Jorgensen  method 
In  this  report. 

The  Newtonian  correction  factor  Is  only  valid  in  the  Hypersonic  Mach  re¬ 
gime.  Some  limited  empirical  results,  ratioed  to  Newtonian  values,  are  avail¬ 
able  and  plotted  in  Figure  34.  This  ratio,  K  ,  is  Included  as  a  multi- 

n 

plier  to  the  viscous  component  of  both  the  normal  force  and  pitching  moment 
coefficients.  The  accuracy  of  this  correction  factor  is  demonstrated  in 
Figure  35.  The  dashed  lines  represent  the  use  of  the  Integral  form  of  the 
equations,  given  above;  the  solid  lines  represent  the  same  predictions  with 
the  correction  term  added  to  correct  the  Newtonian  theory.  The  latter  method 
Is  demonstrated  to  be  superior. 

3. 2. 1.3  Lateral -01 rectional  Coefficients  -  Since  the  Jorgensen  method  is  used 

for  predicting  Ckl  and  C  ,  the  same  approach  was  selected  for  and 
N  m  Y 

C  .  The  slender  body  and  Newtonian  theory  correction  terms  become  roll 
dependent.  Hence,  the  selection  of  the  proper  roll  attitude  and  axis  system 
is  required  to  obtain  the  correct  coefficients. 

3.2.2  Arbitrary  Bodies  with  Nonelllptlc  Cross  Sections 

The  capability  to  predict  the  linear  portion  of  the  longitudinal  stability 
coefficients  of  bodies  with  arbitrary  cross  sections  has  been  included  in  the 
code.  The  methods  selected  for  predicting  normal  force  and  pitching  moment 
curve  slopes  are  based  on  slender  body  theory  since  it  allows  for  rapid  pre¬ 
dictions  without  requiring  the  user  to  panel  the  configuration.  Slender  body 
theory  also  has  the  versatility  of  being  Mach  number  Independent.  To  predict 
the  pressure  drag  of  arbitrary  bodies,  a  panel  model  of  the  body  is  created 
from  the  user  inputs.  The  code  employed  is  based  on  the  Supersonic/Hypersonic 
Arbitrary  Body  Program  (S/HABP). 
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EXPERIMENT 


Even  though  the  approach  selected  did  allow  for  the  analysis  of  a  wide 
array  of  body  cross  sections,  the  basis  In  slender  body  theory  resulted  In 
some  limitations  on  the  types  of  geometries  which  could  be  accurately  ana¬ 
lyzed.  It  was  found  that  the  nose  flness  ratio  had  to  be  greater  than  five. 
The  body  base  cross-sectional  aspect  ratio  had  to  be  less  than  five.  The  base 
cross-sectional  apect  ratio  was  defined  as  the  base  width  squared  and  divided 
by  the  base  area.  The  final  restriction  Is  that  the  bodies  can  not  be  boat- 
tailed  or  flared. 

The  following  sections  describe  the  methods  used  to  predict  the  indivi¬ 
dual  coefficients.  Currently  methods  are  available  to  predict  CM,  C  , 

N  m 

C„,  C  and  C.  .  Notably  no  simple  methods  such  as  those  described  In 

Y  n  Aa 

previous  sections  have  been  found  for  estimating  the  values  of  Ckl  ,  C„  , 

No  No 

Cv  ,  CA  . 

Yo  Aa 

3.2.2. 1  Axial  Force  -  The  axial  force  methods  implemented  In  the  code  vary 

with  the  free  stream  Mach  numbers  being  evaluated.  There  Is  one  comb1r~tion 

of  methods  for  the  subsonic/transonic  regime  and  another  comblnatl  of 

methods  for  the  supersonic  regime.  The  C.  predictions  are  limited  to  zero 

A 

degrees  angle  of  attack  and  consist  of  a  skin  friction  and  a  pressure  drag 
component: 


The  term  Is  the  skin  friction  component  which  Is  proportional  to  the 
wetted  area.  The  pressure  term  is  proportional  to  the  body  fineness  ratio. 

3. 2. 2. 1.1  Skin  Friction  -  The  skin  friction  portion  of  the  axial  force  co¬ 
efficient  Is  computed  using  the  circular  body  skin  friction  methods.  These 
methods  compute  the  flat  plate  skin  friction  coefficient  for  the  free  stream 
Reynolds  number.  The  flat  plate  skin  friction  coefficient  is  then  corrected 
for  curvature  and  non-dimensional i zed  by  the  ratio  of  the  wetted  area  divided 
by  the  reference  area. 
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Quadra laterals  are  computed  using  the  S/HABP  geometry  package  to  approxi¬ 
mate  the  surface.  The  wetted  area  Is  the  sum  of  the  area  of  the  quadralater- 
als. 

3. 2. 2. 1.2  Pressure/Wave  Drag  -  At  a  free  stream  Mach  numbers  of  less  than 
1.4,  a  subsonic/transonic  pressure  drag  method  is  used.  This  method  uses  the 
simple  expression  for  three-dimensional  pressure  drag  presented  in  the  feas¬ 
ibility  study. 

CA  *  7(t/C)3  cf  Sc 
p  ^Ref 

-  body's  equivalent  axlsymmetrlc  diameter 

-  body's  length 
-■  skin  friction  coefficient 

-  cross  sectional  area 

-  reference  area 

v 

The  accuracy  of  this  method  has  not  been  checked  for  a  wide  range  of  con¬ 
figurations.  It  was  developed  for  ellipsoid  type  bodies  which  as  a  class  of 
shapes  which  have  low  cross  sectional  aspect  ratios.  The  small  number  of  com¬ 
parisons  is  due  to  a  lack  of  subsonic  arbitrary  body  data.  The  lack  of  data 
may,  however,  reflect  a  lack  of  interest  in  arbitrarily  shaped  subsonic  con¬ 
figurations.  Further  method  validation  and  development  is  necessary. 

For  free  stream  Mach  numbers  greater  than  1.4,  the  pressure  drag  is  com¬ 
puted  using  S/HABP  methods.  The  nose,  which  Is  described  by  the  user  using 
the  nose  tip  and  base  cross  section  shapes.  Is  broken  up  Into  quadralaterals . 
The  pressure  coefficients  on  each  of  these  panels  is  computed  using  the  ACM 
empirical  equations  for  leeside  panels  and  the  Dahlem-Buck  equations  for  wind¬ 
ward  panels  (Reference  18).  The  pressure  force  on  each  panel  Is  then  com¬ 
puted  and  resolved  into  the  axial  force  direction.  The  resulting  axial  force 
is  then  nondimensionalized  to  compute  the  portion  of  the  axial  force  coeffi¬ 
cient  which  results  from  the  pressure  forces. 


59 


The  pressure  forces  are  computed  only  for  the  nose  in  the  arbitrary  body 
methods.  This  Is  because  the  body  is  assumed  to  be  at  zero  degrees  angle  of 
attack.  Since  the  centerbody  sides  are  parallel  to  the  centerline,  it  will 
not  contribute  pressure  forces  in  the  axial  direction. 

Comparison  of  the  axial  force  coefficient  results  obtained  using  this 
method  with  data  showed  an  overprediction  trend.  The  median  of  the  predicted 
values  was  found  to  be  25  percent  above  the  data  (Figure  36).  Therefore,  a 
correction  factor  was  applied  to  reduce  the  predicted  axial  force  oy  25  per¬ 
cent. 


3. 2. 2. 2  Normal  Force  Coefficient  -  The  methods  coded  to  predict  the  normal 

force  coefficients  (C.,)  of  arbitrary  cross  section  shaped  bodies  are  based 

to 

on  slender  body  theory.  Methods  have  been  coded  to  predict  only  the  normal 
force  coefficient  slope.  The  normal  force  at  zero-angle-of -attack  Is  not 
currently  computed.  Currently  no  methods  are  included  to  predict  the  viscous 
contribution  to  normal  force  coefficient. 


The  accuracy  of  the  selected  method  is  based  on  tuw  ability  to  accurately 
evaluate  the  apparent  area  of  the  cross  section  at  the  base  of  the  nose.  A 
modified  method  of  Hess  and  Sm  th,  (Reference  19)  was  used  to  calculate  appar¬ 
ent  area.  The  nose  cross  section  Is  divided  into  100  segments,  each  with  a 
distributed  source  strength.  The  results  of  this  method  have  been  checked 
against  analytical  solutions.  With  100  test  points  the  numerical  results  were 
within  three  percent  of  the  analytical  results  for  all  cross  sections  except 
thin  vertical  shapes.  Figure  37  shows  comparisons  of  the  method  used  to  the 
analytical  closed  form  solutions  for  rounded  shapes,  and  Figure  38  shews  the 
same  comparisons  for  flat  sided  shapes. 


Once  the  apparent  area  of  the  cross  section  has  been  computed,  the  normal 
force  coefficient  is  computed  from  the  equation: 


where: 


CNa  =  2JV 
$Ref 

-  apparent  area  of  base  cross  section  of  the  nose 

-  reference  area 
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FIGURE  37  COMPARISON  OF  ADDED  MASS  VALUES  FOR.  VARIOUS  ROUNDED  CROSS  SECTIONS 
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FIGURE  38  COMPARISON  OF  ADDED  MASS  VALUES  FOR  VARIOUS  FLAT  SIDED  CROSS  SECTIONS 


The  simplicity  of  this  equation  means  that  the  value  of  normal  force  co¬ 
efficient  Is  Mach  number  Independent  and  only  a  function  of  the  nose  base 
cross  section.  However,  there  are  some  limitations  placed  on  the  configura¬ 
tions  for  which  accurate  results  can  be  obtained  due  to  the  methods  basis  In 
slender  body  theory.  Comparison  with  data  have  shown  that  the  fineness  ratio 
of  the  nose  should  be  greater  than  five.  The  aspect  ratio  of  the  cross  sec¬ 
tion: 


«c.s- 


WIDTH  2 

mnm 


must  be  less  than  five.  These  two  llmlations  keep  the  configuration  from  be¬ 
coming  too  blunt  or  flat. 


To  calculate  the  contribution  of  the  centerbody  which  Includes  the  pres¬ 
sure  carryover  from  the  nose  an  equivalent  circular  body  method  has  been  em¬ 
ployed.  An  equivalent  axisymmetrlc  body  radius  Is  computed  based  on  the  maxi¬ 
mum  arbitrary  body  cross  section.  The  relationship  used  Is: 


Using  the  value  of  r  ,  the  C..  of  a  circular  cone-cylinder  conflgura- 

eq  Na 

tlon  with  the  same  nose  and  centerbody  fineness  ration  as  the  arbitrary  body 

is  computed  using  the  Missile  Datcom  axisymmetrlc  body  methods.  This  process 

is  outlined  In  Figure  39  for  a  rounded  triangular  body.  After  the  equivalent 

cone-cylinder  C  is  computed.  It  is  scaled  based  on  the  ratio  of  the  ap- 
Na 

parent  areas  of  the  arbitrary  body  to  that  of  a  cone.  The  equation  that  re¬ 
sults  Is: 

CNoARB.B0DY  *  A'  ARB, BODY  %,  rQNE/CYL 
A'  CONE 


The  results  of  this  method  have  been  checked  against  data  and  the  results 

of  a  comparison  with  one  extensive  data  base  are  presented  in  Figure  40.  The 

configurations  In  this  data  base  were  tested  at  five  Mach  numbers  from  2.5  to 

4.5.  Figure  40  presents  the  error  in  the  prediction  of  C  for  each 

N 
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EXP.  CALC  ERROR 
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FIGURE  40  COMPARISON  OF  Cn  DATA  FOR  VARIOUS  NOSE/CENTERBODY 
CONFIGURATIONS  OF  FINENESS  RATI0=10 


configuration  tested  at  a  Hach  number  of  2.5.  The  worst  case  is  the  2.5  to  1 
horizontal  ellipse  with  an  error  of  16.2  percent.  From  this  figure  it  can  be 
seen  that  the  accuracy  is  acceptable  for  preliminary  design. 

This  method  Is  considered  conducive  to  preliminary  design  since  it  is  both 
quick  running  and  allows  for  the  body  cross  sections  to  be  changed  easily. 
These  features  allow  for  the  types  of  parametric  studies  commonly  performed  in 
the  preliminary  design  environment. 

3. 2. 2. 3  Pitching  Moment  Coefficient  -  As  with  the  normal  force,  only  the 
linear  component  of  the  pitching  moment  is  computed.  There  are  no  methods 
are  Included  for  computing  either  zero  angle  of  attack  pitching  moment  or  the 
nonlinear  component. 

The  method  incorporated  Into  the  code  for  computing  C  is  based  on 

the  apparent  area,  A1,  calculation  used  to  compute  .  Using  slender 

body  theory  the  value  of  C  is  found  from  the  Integral  equation. 

me* 


base 

=  /  _ ? _  .  X  •  dA' 

J  sRofLRef 
nose 


were  C  Is  measured  about  the  nose  tip. 

Ha 

This  equation  is  automated  by  breaking  the  nose  up  into  nine  longitudinal 
stations  for  which  the  apparent  corss-sect.ional  areas  are  computed.  The 
C  contribution  of  each  segment  is  computed  using  the  equation: 

moL 


u.  " - (*cg  “  xi)  ( '  -  A-j-i ) 

1  ^Refl-Ref 
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■ratted- 


where: 


A  *  1 

A'i-i 

LRef 

SRef 

XCE 


-  apparent  area  of  segments  base  cross  section 

-  apparent  area  of  segments  front  cross  section 

-  reference  length  (usually  base  diameter) 

-  reference  area  (usually  base  area) 

-  longitudinal  center  of  gravity  position 

-  centroid  of  1th  segment 


The  contributions  of  each  segment  to  the  nose  C  Is  then  summed. 

m& 


The  pitching  moment  curve  slope  value  that  is  computed  for  the  nose  is 
corrected  to  Include  the  effect  of  an  afterbody  by  using  a  method  analogous 
to  the  normal  force  curve  slope  correction.  The  method  employed  uses  the  com¬ 
putation  of  an  equivalent  axlsymmetrlc  body  to  provide  the  nose  to 

afterbody  ratio  of  C  effects.  Oata  comparisons  of  nose  alone  results 

m& 

have  shown  sufficient  accuracy  for  preliminary  design  analysis. 
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4.0  FIN  ALONE  METHOD  SELECTIONS 


The  method  summary  for  fin  alone  aerodynamics  is  shown  in  Figure  41.  Un¬ 
like  the  method  selected  for  body  alone  aerodynamics,  most  of  the  subsonic/ 
transonic  methods  are  theoretical  or  vemi empirical .  These  types  of  methods 
maximize  the  panel  shapes  that  can  be  addressed. 

Table  2  summarizes  the  methods  selected  compared  to  those  recommended  In 
the  Feasibility  Study.  Note  that  most  of  the  methods  implemented  are  those 
recommended  during  the  Feasibility  Study.  The  three  methods  changed  were  due 
to  ease  of  use,  extended  capability,  or  additional  accuracy  through  use  of  a 
theoretical  method. 

Figure  42  Indicates  those  fin  method  application  areas  as  a  function  of 
Mach  regime  and  panel  shape.  Note  that  the  area  with  the  most  questionable  or 
limited  application  Is  that  of  ion -straight  tapered  fins.  Although  airfoil 
section  results  are  not  generally  available  at  supersonic  speeds,  this  Is  not 
a  major  problem  area  since  the  methods  at  these  speeds  assumed  that  the  fins 
are  thin.  Supersonic  fins  usually  are  thin  to  minimize  drag. 

The  following  paragraphs  describe  the  method  selections  for  straight 
tapered  and  nonstraight  tapered  fin  panels. 

4.1  FINS  WITH  STRAIGHT  TAPER 


There  are  far  more  methods  for  predicting  the  aerodynamics  of  straight 
tapered  fins  than  there  are  for  other  panel  shapes.  This  is  particularly  true 


at  transonic  and  supersonic  speeds, 
methods  for  fins  than  for  any  other  component  of  the  configuration,  and  this 
Is  probably  due  to  the  amount,  of  work  performed  to  analyze  airplanes,  parti¬ 
cularly  commercial  aircraft. 


ill 
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FIGURE  41  FIN  ALONE  METHOD  SUMMARY 
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The  following  paragraphs  describe  the  methods  selected  for  straight  taper¬ 
ed  fins,  typical  of  those  used  on  conventional  missile  configurations.  Note 
that  methods  were  selected  which  were  applicable  to  panels  with  aspect  ratios 
from  0.5  to  4.  Although  panel  aspect  ratios  greater  than  4  have  been  predict¬ 
ed  with  the  selected  methods  yielding  excellent  results,  they  have  not  been 
verified  for  a  wide  variety  of  panels.  Therefore,  these  methods  can  not  be 
used  with  a  high  degree  of  confidence  for  high  aspect  ratio  panels. 

4.1.1  Axial  Force  Coefficient 

As  with  bodies  alone,  the  axial  force  of  a  fin  Is  assumed  to  be  comprised 
of  a  zero-lift  component  and  a  component  due  to  angle  of  attack.  However  fins 
can  be  cambered,  which  can  result  In  the  zero-lift  component  of  axial  force 
not  occurlng  at  zero  angle  of  attack.  Therefore,  the  methods  for  fin  alone 
aerodynamics  are  more  complex,  and  often  require  a  method  change  as  a  function 
of  Nach  number  and  angle  of  attack. 

Figure  43  Illustrates  the  variation  In  fin  drag  (or  axial  force)  at  zero- 
lift.  The  subsonic  pressure  drag  method  assumes  that  It  Is  almost  Invariant 
to  a  Nach  number  of  1.05,  and  then  1  Inearl ly  decreases  to  zero  at  Nach  1.2 
Wave  drag  Is  assumed  to  begin  at  Nach  0.8,  when  local  shocks  first  appear  In 
the  flow,  Increasing  to  a  maximum  near  Nach  1.Q5,  where  the  flow  Is  all 
supersonic.  Leading  edge  bluntness  drag  varies  nearly  monotonlcally  across 
the  Nach  spectrum. 

4. 1.1.1  Skin  Friction  -  The  skin  friction  of  fins  can  be  calculated  using  the 
method  described  for  body  alone,  except  that  a  two-dimensional  to  three-dimen¬ 
sional  correction  factor  to  account  for  surface  shape  is  not  used.  This  is 
due  to  the  fact  that  the  skin  friction  methods  are  based  upon  flat  plate  data 
correlations  and  theory.  Although  experimental  data  for  skin  friction  is  dif¬ 
ficult  to  obtain,  the  accuracy  of  the  method  used  has  been  Implied  from  numer¬ 
ous  correlations  with  axial  force  data  across  the  Nach  regime. 
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43  FIN  DRAG  METHODS 


4. 1.1. 2  Subsonic  Pressure  Drag  -  The  subsonic  assumption  implies  that  no 
shocks  are  present  In  the  flow  field  region  being  analyzed.  In  the  transonic 
flow,  a  particular  fin  panel  can  have  both  subsonic  and  supersonic  flow 
occurring  simultaneously  over  different  portions  of  the  panel.  For  bodies,  it 
was  assumed  the  subsonic  pressure  drag  was  a  function  of  the  body  length  to 
diameter  ratio,  and  was  correlated  to  the  configuration  skin  friction  coeffi¬ 
cient.  This  scheme  is  used  for  fins  as  well,  except  that  the  correlation  is 
not  based  upon  overall  length,  but  rather  the  thickness  to  chord  ratio. 

4. 1.1. 3  Transonic  Pressure/Wave  Drag  -  This  nomenclature  is  used  for  that 
portion  of  fin  drag  which  exists  due  to  presence  of  shocklets  In  the  flow 
field.  This  occurs  near  Mach  0.8  and  Increases  until  the  fin  flow  Is  entirely 
supersonic  (except  perhaps  blunted  leading  edges  where  shocklets  may  be  pre¬ 
sent  at  lower  Mach  numbers).  The  method  used  is  a  potential  flow  method, 
(Reference  6),  and  is  solved  by  calculating  the  potential  flow  about  the  fin 
at  Mach  1.05,  and  then  llnearlly  reducing  that  value  as  a  function  of  Mach 
number  to  zero  at  Mach  0.8.  Although  this  approach  Is  approximate,  the  re¬ 
sults  have  been  observed  to  give  accuracies  consistent  with  preliminary  design. 

4. 1.1. 4  Supersonic  Wave  Prag  -  The  methods  available  for  calculating  fin 
supersonic  wave  drag  have  been  developed  for  thin  airfoil  sections.  Correc¬ 
tions  are  available  to  extrapolate  the  method  to  surfaces  with  small  to  mode¬ 
rate  thickness  to  chord  ratios.  The  extrapolation  Is  developed  primarily 
from  experimental  observations.  The  method  selected  was  developed  for  the 
NSWC  Aeropredlctlon  code  (Reference  29),  where  the  surface  of  the  fin  Is 
transformed  into  a  computational  grid,  and  a  potential  flow  solution  is 
applied.  This  method  is  very  powerful,  and  permits  the  accurate  calculation 
of  the  drag  characteristics  of  straight-tapered  surfaces,  but  has  the  dis¬ 
advantage  that  no  source  or  sink  line  placed  on  the  airfoil  surface  can  be 
swept  forward;  this  condition  is  easily  violated  if  the  surface  has  a  forward 
swept  trailing  edge.  Hence,  both  this  method,  and  that  f  ~om  the  U.  S.  Air 
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Force  DATCOM  (Reference  30)  which  allows  for  swept  forward  trailing  edges, 
are  Incorporated  Into  the  computer  code.  The  potential  flow  method  is  more 
accurate  than  the  DATCOM  method. 

4. 1.1. 5  Leading  Edge  Bluntness  Drag  -  The  supersonic  source/sink  method 
described  In  Section  4. 1.1. 4  has  the  capability  to  compute  the  effect  of 
spherically  blunted  leading  and  trailing  edges  on  fin  panels.  At  subsonic  and 
transonic  speeds,  however,  there  Is  no  easy-to-apply  accurate  method  avail¬ 
able.  Hence,  the  empirical  correlation  from  the  U.S.  Air  Force  DATCOM,  which 
Is  a  function  of  leading  edge  projected  area,  Mach  number,  and  leading  edge 
sweepback  was  selected.  There  was  Insufficient  data  to  completely  evaluate 
the  accuracy  of  the  technique,  but  it  is  believed  to  be  sufficiently  accurate 
for  missile  preliminary  design. 

4. 1.1. 6  Base  Drag  -  The  method  selected  for  fin  base  drag  Is  based  on  two- 
dimensional  base  drag  data.  This  method  could  not  be  totally  validated  since 
the  base  drag  of  fins  is  normally  such  a  small  contribution  to  total  drag. 
However,  the  data  used  in  deriving  the  method  is  believed  to  be  accurate  and 
since  this  component  is  small,  It  does  not  warrent  a  high  degree  of  validation. 

4. 1.1. 7  Angle-of-Attack  Effect  -  As  mentioned  earlier,  the  drag  of  a  fin  at 
angle  of  attack  has  been  extensively  studied,  both  theoretically  and  empiri¬ 
cally,  for  aircraft  configurations.  The  method  selected  for  Missile  DATCOM  Is 
that  of  the  U.  S.  Air  Force  DATCOM.  The  plot  of  drag  versus  lift  Is  assumed 
to  have  a  parabolic  shape  and  to  be  dependent  of  fin  panel  lift,  airfoil  sec¬ 
tion  shape,  and  aspect  ratio  efficiency  (Oswald's  efficiency).  Results  com¬ 
puted  using  this  method  have  agreed  quite  well  with  experiment  for  a  wide 
variety  of  missile  configurations.  However,  the  method  seems  to  fall  when  the 
angle  of  attack  approaches  30  degrees;  the  drag  due  to  lift  curve  shape  de¬ 
parts  significantly  from  the  parabolic  shape.  This  method  discrepancy  becomes 
more  noticeable  for  missiles,  where  the  required  angle-of-attack  range  is  much 
higher  than  for  airplanes.  It  is  recommended  that  a  new  method  be  developed 
which  addresses  the  higher  angle-of-attack  requirements  of  missiles. 
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4.1.2  Normal  Force  and  Pitching  Moment  Coefficient 


Figure  44  Illustrates  the  method  used  to  compute  fin  normal  force  as  a 
function  of  angle  of  attack  and  Mach  number.  The  normal  force  Is,  analagous 
to  that  for  bodies  alone.  It  Is  assumed  to  be  comprised  of  a  potential  and  a 
nonlinear  contribution: 


CN  “  CNP  +  CNN 

where, 

CNP  *  CNaSln<"> 


and 


'NN 


W"  <“> 


The  potential  normal  force  is  computed  from  the  fin  panel  C  ,  which 

Na 

is  Mach  dependent.  The  nonlinear  effect  Is  analagous  to  the  viscous  effect 
for  bogles,  except  that  It  is  assumed  to  be  a  stall  or  post  stall  effect.  An 
alternate  method  considered  uses  a  "true"  viscous  model  where  the  effect  of 
angle-of -attack  Is  assumed  to  be  is  obtained  at  90  degrees  angle  of  attack, 
adjusting  for  cross  flow  Mach  number.  This  Idea  was  not  Implemented  since  it 
Ignores  the  effect  of  stall. 


A  key  Input  Into  the  angle-of -attack  analysis  is  the  angle  of  attack  for 
stall  aM  shock  detachment.  As  shown  In  Figure  44,  the  method  used  is  depen¬ 
dent  on  Mach  number,  angle-of-attack,  the  airfoil's  thickness  to  chord  ratio, 
and  lift. 


4. 1.2.1  Airfoil  Section  Charactarlst'cs  -  At  subsonic  and  transonic  speeds, 
the  airfoil  section  characteristics  are  required  in  order  to  accurately  model 
fin  panel  lift  and  drag.  The  method  selected  for  incorporation  Into  the  com¬ 
puter  code  Is  the  Airfoil  Section  Module  from  the  Digital  DATCOM  computer 
code,  which  Is  used  for  the  prediction  of  aircraft  aerodynamics.  It  has  been 
validated  for  a  wide  variety  of  airfoil  shapes  and  is  a  relatively  simple 
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technique  to  apply.  However,  It  Is  not  amenable  to  hand  calculation,  and  the 
Handbook  user  must  rely  upon  data  from  an  airfoil  section  handbook  to  obtain 
the  airfoil's  aerodynamic  characteristics.  The  airfoil  section  characteris¬ 
tics  of  the  fins  are  not  required  at  supersonic  speeds,  since  those  methods 
assume  thin  airfoils,  and  a  limited  correction  for  thickness  to  chord  ratio  is 
used. 

4. 1.2. 2  Potential  Normal  Force  and  Hitching  Moment  Coefficient  -  Since  the 
potential  pitching  moment  coefficient  Is  computed  using  the  potential  normal 
force  coefficient.  It  Is  important  that  the  potential  normal  force  be  as  accu¬ 
rate  as  possible.  In  order  to  Insure  accuracy  over  the  greatest  range  of 
planform  shapes,  results  from  theoretical  methods  were  selected  If  possible  at 
all  Mach  numbers.  At  subsonic  speeds,  the  most  reliable  method  is  the  Lowry- 
Polhamus  correlation  for  straight  tapered  surfaces  (Reference  31).  It  Is  a 
theoretical  method  which  had  been  adjusted  to  correlate  with  experimental 
data.  It  has  been  validated  over  wide  range  of  panel  aspect  ratios,  down  to 
an  aspect  ratio  of  one  half.  Its  accuracy  Is  unknown  for  aspect  ratios 
smaller  than  0.5;  no  upper  limit  has  been  noted  in  the  method  validation. 

At  transonic  speeds,  the  only  theoretical  results  available  are  those 
from  the  R.A.S.  data  sheets  (Reference  32).  The  curves  from  the  R.A.S.  com¬ 
pendium  represent  results  for  panels  with  aspect  ratios  that  cannot  be  accu¬ 
rately  evaluated  in  the  transonic  regime,  but  are  "best"  estimates  based  upon 
experimental  data  results.  To  use  a  true  theoretical  method  at  transonic 
speeds  would  be  cost  prohibitive,  hence  an  Interpolation  of  this  data  base  is 
the  best  compromise  of  cost  and  accuracy.  The  predictions  coincide  with  the¬ 
oretical  results  at  both  subsonic  and  supersonic  speeds  for  method  compat¬ 
ibility. 

At  supersonic  speeds,  supersonic  wing  theory  was  selected.  The  results 
have  been  shown  to  be  accurate  for  a  wide  range  of  thin  wings  which  makes  the 
method  a  powerful  missile  aerodynamic  prediction  tool. 
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4. 1.2. 3  Nonlinear  Normal  Force  and  Pitching  Moment  Coefficient  -  The  non¬ 
linear  variation  In  normal  force  Is  based  upon  the  method  of  the  U.S.  Air 
Force  DATCOM.  The  method  Is  an  empirical  correlation  for  wing  normal  force 
coefficient  with  angle -of -attack  for  a  wide  variety  of  wing  shapes.  Although 
the  form  of  the  method  Is  the  same  at  all  speeds,  the  phenomena  modeled  are 
quite  different.  At  subsonic  speeds,  the  wing  stall  angle  of  attack  and  wing 
maximum  lift  coefficient  (at  stall)  are  required;  these  are  obtained  using 
the  U.S.  Air  Force  DATCOM  methods  as  well.  The  stall  characteristics  dominate 
the  post  stall  behavior  of  the  wing  lift  coefficient.  Since  this  post  stall 
behavior  is  based  upon  empirical  data  correlations  for  "typical1'  subsonic 
airfoils,  the  characteristic  shape  of  the  wing  stall  will  not  necessarily  be 
that  experienced  for  thin,  missile  airfoils.  Since  no  other  method  Is 
readily  available,  this  method  was  selected.  Beyond  stall,  the  wing  lift  for 
missile  wings  closely  follows  results  typical  of  aircraft. 

At  supersonic  speeds,  the  concept  of  stall  Is  not  directly  applicable. 
Instead,  the  wing  lift  will  vary  almost  llnearlly  until  the  angle-of-attack  Is 
great  enough  for  the  leading  edge  shock  to  become  detached.  The  shock  detach¬ 
ment  angle  of  attack  Is  a  function  of  the  Hach  number  and  the  semi-apex  angle 
of  the  airfoil  leading  edge.  Above  shock  detachment,  the  U.S.  Air  Force  DAT¬ 
COM  methods  for  wing  lift  variation  Is  used.  The  shock  detachment  angle  Is 
determined  using  the  DATCOM  method. 

At  transonic  speeds  there  Is  no  suitable  method  for  determining  the  angle 
of  attack  for  which  wing  lift  no  longer  remains  linear.  An  approximate  method 
was  Implemented.  It  assumes  that  the  maximum  lift  angle  of  attack  In  shock- 
free  flow  ( Subsonic  soeedsl  and  +  ho  chnrk  rfptarhn»nt  annlp  nf  attjirk  (at 

»  r  C  ~  ~  *•  O  *  ~  “  ■  - -  ^  —  - 

supersonic  speeds)  can  be  llnearlly  Interpolated  as  a  function  of  Hach  number 
at  transonic  speeds.  Results  obtained  using  this  procedure  have  been  found  to 
have  adequate  accuracy  for  preliminary  design. 

Figure  45  presents  the  range  of  planform  shapes  and  Mach  numbers  used  In 
the  validation  of  fin  alone  normal  force  methods.  The  aspect  ratios  Investi¬ 
gated  where  from  0.5  to  4.0,  the  taper  ratios  from  0.0  to  1.0,  and  the  Hach 
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FOR  WHICH  FIN  ALONE  ~N  PREDICTIONS  WERE  EVALUATED 


numbers  from  0.80  to  4.60.  Typical  results  are  shown  In  Figure  46  for  an  as¬ 
pect  ratio  0.5  panel  with  a  nonswept  trailing  edge.  The  four  Mach  numbers 
shown  Include  a  subsonic/transonic,  a  low  supersonic,  and  two  supersonic  con¬ 
ditions.  Up  to  30°  angle  of  attack,  the  prediction  method  follows  the  ex¬ 
perimental  data  closely.  Similar  results  were  obtained  for  all  of  the  fin 
panels  Illustrated  In  Figure  45. 


The  potential  pitching  moment  Is  computed  using  the  linear  center  of  pres¬ 
sure,  which  was  obtained  using  theoretical  methods  similar  to  those  used  for 
normal  force  coefficient.  The  nonlinear  center  of  pressure  was  assumed  to  be 
at  the  panel  area  centroid  at  90°  angle  of  attack  for  all  speeds.  The 
center  of  pressure  was  assumed  to  vary  linearily  as  a  function  of  angle-of- 
aitack  from  the  linear  center  of  pressure  at  zero  angle  of  attack  to  the  area 


81 


rOIICS  COEFFICIENT  NORMAL  FORCE  COEFFICIENT 


11-4*01 

lo-aa-sa 


FIGURE  46  FIN  ALONE  Cn  OATA  COMPARISON  FOR  AN  ASPECT  RATIO  0.5  FIN 


82 


centroid  at  90°  angle  of  attack.  This  method  is  an  engineering  approxima¬ 

tion  since  no  theoretical  method  which  can  be  easily  applied  is  available. 

4.1.3  Lateral -Directional  Coefficients 

Specific  methods  for  computing  the  lateral-directional  characteristics  of  fin 
panels  alone  were  not  selected  for  application  in  Missile  DATCQH.  The 
lateral-directional  coefficients  are  computed  for  complete  configurations 
using  a  concept  called  the  "equivalent  angle-of -attack"  method  which  uses  the 
panel  normal  force  and  center  of  pressure  normal  to  the  panel.  This  Is  ad¬ 
justed  to  account  for  panel  deflection  and  roll  angle  about  the  fuselage. 

This  method  yields  the  panel  loads  normal  to  the  panels  which  are  then  re¬ 

solved  Into  body  axis  components.  Therefore,  in  the  calculation  of  configura¬ 
tion  normal  force  and  pitching  moment,  the  loads  for  computing  side  force, 

yawing  moment,  and  rolling  moment  are  directly  available.  This  method,  does 
not  require  another  method  dedicated  to  the  calculation  of  lateral-directional 
coefficients. 

4.2  FINS  WITH  NON-STRAIGHT  TAPER 

The  methods  for  non-straight  tapered  fins  are  essentially  those  for  straight 
tapered  fins.  The  exceptions  are  the  calcuatlon  of  normal  force  coefficient 
and  aerodynamic  center. 

For  normal  force  and  normal -force-curve  slope,  an  effective  straight 
tapered  panel  is  computed.  This  effective  panel  has  the  same  aspect  ratio 
and  taper  ratio  as  the  nor.straight  panel  but  uses  an  effective  h  "f-chord 
sweep  angle.  This  angle  Is  computed  as: 

n 

coslAc/2)eff  *  r  £  cosUc/2)i  SF. 

i  =  l 
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where 


N  Is  the  number  of  panel  segments 

Sp  is  the  total  panel  area 

(Ac/2),.  is  the  half-chord  sweep  angle  of  each  panel  segment 

5r,  Is  the  area  of  each  panel  segment 

r  1 

Normal  force  and  normal -force-curve  slope  coefficients  are  then  computed  using 
this  effective  straight-tapered  panel. 

To  compute  the  panel  aerodynamic  center  for  nonstraight  tapered  fins,  the 
fin  Is  divided  Into  two  fin  panels.  The  individual  lift  and  aerodynamic 
center  for  each  panel  are  then  used  to  establish  the  aerodynamic  center  for 
the  complete  fins. 
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5.0  INLET  METHODS 


The  inlet  methods  incorporated  into  Missile  DATCOM  are  designed  to  provide 
the  capability  to  analyze  the  external  aerodynamics  of  general  Inlet  geome¬ 
tries  without  the  need  for  complex  inputs.  This  capability  had  been  developed 
In  a  pilot  program  written  by  Krleger,  Williams  and  Hood  (Reference  37). 

The  pilot  program  was  restricted  to  the  analysis  of  axisymmetric  and  2-D 
inlets  operating  on  vehicles  flying  at  supersonic  speeds.  Only  the  external 
forces  on  the  Inlet  are  predicted  and  the  assumption  is  made  that  the  Inlet  is 
operating  at  its  design  mass  flow.  No  spill  drag  or  body  inlet  interference 
effects  were  Included  in  the  methods. 

This  pilot  program  was  used  as  the  basis  for  the  Missile  DATCOM  methods. 
No  methods  were  identified  that  provided  accurate  subsonic  or  transonic  pre¬ 
dictions.  As  a  result,  no  capability  to  analyze  Inlets  operating  in  subsonic 
or  transonic  flow  was  incorporated  Into  Missile  DATCOM. 

The  contributions  of  the  inlets  to  the  full  configuration  aerodynamics  are 
treated  as  Incremental  additions.  Up  to  four  axisymmetric  or  2-D  inlets 
mounted  on  diverters  can  be  positioned  around  the  missile  body  at  user-defined 
roll  angles.  The  contribution  of  the  body  surface  covered  by  the  diverter 
is  subtracted  from  the  body  coefficients  In  which  it  had  been  included.  Both 
the  skin  friction  and  pressure  forces  are  accounted  for  on  the  surfaces  of  the 
inlet.  No  other  body-inlet  Interference  effects  are  included. 

5.1  METHODOLOGY  USED 

The  inlet  and  diverter  geometry’s  are  input  separately  in  a  coordinate 
system  attached  to  the  inlet.  This  process  is  shown  as  steps  1  and  2  of 
Figure  47.  Figure  47  describes  the  process  by  which  the  inlet  model  is  con¬ 
structed  and  then  its  position  is  located  relative  to  th»  body. 
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ACTUAL  CONFIGURATION 


FIGURE  47  GEOMETRY  CONSTRUCTION  PROCEDURE 


Once  the  Inlet  and  diverter  geometries  have  been  described,  a  panel  model 
of  them  is  created.  This  model  consists  of  S/HABP  quadralaterals  constructed 
by  geometry  package  described  in  Reference  37.  The  corner  points  of  the  quad¬ 
rilaterals  are  defined  using  the  Inlet  diverter  inputs.  A  description  of 
these  Inputs  and  assumptions  made  by  the  code  are  Included  in  Volume  1  of  the 
Mislle  DATCOM  User's  Manual.  The  paneled  geometries  are  shown  in  Figure  48. 
The  axlsymmetric  inlets  are  paneled  circumferentially  while  the  2-D  inlets  are 
constructed  of  longitudinal  panels.  Both  are  attached  to  a  diverter  which  is 
attached  to  the  missile  body. 

The  coefficient  contributions  resulting  from  the  surfaces  of  the  diverter 
that  attach  to  either  the  inlet  or  the  missile  body  are  subtracted  out  of  the 
aerodynamic  calculations  when  the  inlet  increments  are  computed.  This  is  done 
to  negate  the  contributions  computed  for  the  hidden  surfaces  of  the  Inlet  and 
body. 

A  single  Inlet  Is  constructed  In  the  horizontal  position  then  transformed 
into  the  body  axis  system.  This  process  Includes  the  rotation  of  the  Inlet 
about  the  body  centerline  and  the  reflection  of  the  inlet  through  the  vertical 
plane  if  twin  inlets  have  been  prescribed.  Special  consideration  Is  given  to 
twin  Inlets  that  are  input  at  a  rotated  position  of  90  degrees.  For  this  con¬ 
figuration,  the  two  Inlets  are  placed  on  the  top  and  bottom  of  the  body. 

After  the  quadrilaterals  are  transformed,  their  outward  normals  are  com¬ 
puted.  The  angle  between  the  outward  normals  and  the  free  stream  velocity 
vector  is  determined.  This  Is  used  to  compute  the  Impact  angle  of  the  free 
stream  velocity  on  each  panel.  The  pressure  coefficient  is  then  computed  for 
each  panel  individually  based  on  the  Impact  angle.  For  windward  Impact  angles 
the  Dahlem-Buck  equation  is  used,  while  for  leeward  panels  the  ACM  empirical 
equation  is  used.  These  two  equations  are  described  in  Reference  18  including 
a  plot  of  their  values  with  variation  in  impact  angle.  The  pressure  forces  on 
the  inlets  are  computed  by  multiplying  the  area,  the  pressure  coefficient  and 
the  unit  vector  of  each  quadrilateral  then  adding  the  product  to  obtain  the 
total  aerodynamic  coefficient.  The  panel  moment  arm  about  its  reference 
center  of  gravity  is  Included  in  the  product,  for  pitching  moment  coefficient. 
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The  skin  friction  component  of  the  Inlet  forces  is  accounted  for  by  com¬ 
puting  the  average  skin  friction  coeffirient  for  an  equivalent  flat  plate  hav¬ 
ing  equal  length  and  area  of  the  inlet.  The  method  of  Van  Driest  II  is  used 
for  the  skin  friction  coefficient.  The  resulting  shear  force  is  then  correct¬ 
ed  for  three  dimensional  effects.  Skin  friction  effects  are  only  included  in 
the  axial  force  coefficients. 

5.2  RESULTS  AND  IMPROVEMENTS 

The  inlet  method's  prediction  capabilities  were  checked  against  a  data 
base  which  Includes  four  axl symmetric  Inlet  configurations  and  four  two-dimen¬ 
sional  inlet  configurations.  This  data  base  is  documented  in  References  38 
and  39.  Comparisons  were  made  for  all  eight  configurations  at  flight  Mach 
numbers  of  2.5  and  3.95.  The  inlets  were  designed  to  operate  at  a  Mach  number 
of  3.00.  An  oil  flow  analysis  of  the  inlets,  reported  in  Reference  40,  con¬ 
cluded  that  the  axisymmetric  inlet  was  unstarted  at  Mach  2.5  which  agrees  with 
the  design  operation  conditions. 

The  comparisons  showed  good  agreement  at  both  Mach  numbers  and  for 
all  configurations.  Figure  49  shows  a  representative  data  comparison.  It 
presents  the  comparisons  for  the  twin  axisymmetric  configuration  at  Mach  2.5. 
The  Cfto  values  were  found  to  be  acceptable  for  the  higher  Mach  numbers  and 
most  lower  Mach  numbers,  but  the  change  In  with  angle  of  attack  was 
underpredicted  for  most  cases.  The  pitching  moment  prediction  accuracy  was 
found  to  be  very  dependent  on  Mach  number.  The  higher  Mach  number  comparisons 
were  good  (Table  3).  However,  at  the  lower  Mach  number  of  2.5  the  predictions 
were  very  poor.  Figure  49  shows  a  comparison  with  data  that  falls  to  match 
both  the  magnitude  and  trend  of  the  data. 

The  results  of  these  data  comparisons  indicate  that  the  inlet  methods  are 
acceptable  for  supersonic  Mach  numbers  but  only  marginally  for  low  supersonic 
Mach  numbers.  The  poor  predictions  at  the  low  Mach  numbers  are  explainable  by 
the  fact  that  the  inlet  methods  do  not  include  shocks,  vortices  or  the  effects 
of  the  body  flow  field.  The  poor  pitching  moment  comparisons  combined  with 
the  good  normal  force  predictions  indicated  a  local  pressure  distribution 
problem.  A  need  for  Improved  inlet  methods  is  evident. 
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MISSILE  DATCON  CN  PREDICTION  ACCURACY 


MISSILE  DATCOM  INLET  PREDICTIONS  COMPARED  TO  NASA  TN  84558  DATA 


COMPARISON  OF  STATIC  COEFFICIENTS  FOR  BODY  WITH 
2-0  INLET  AT  MACH  =3.95 
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6.0  CONFIGURATION  SYNTHESIS 


The  purpose  of  configuration  synthesis  is  to  take  the  body  and  fin  alone 
aerodynamic  characteristics  and  combine  them  to  obtain  the  total  configuration 
aerodynamic  characteristics.  There  are  many  ways  to  synthesize  a  configura¬ 
tion.  The  most  common  approach  Is  to  take  the  panel  alone  normal  force  co¬ 
efficient,  multiply  It  by  a  carryover  Interference  factor,  obtained  using 
slender  body  theory,  and  then  add  the  result  to  the  body  alone  normal  force. 
Tnls  can  be  represented  mathematically,  for  a  fin-body  configuration,  as: 

CNFB  "  CNB  +  tKF(B)+KB(F)]CNF 

The  disadvantage  of  this  method  Is  that  it  assumes  that  the  angle  of 
attack  range  to  be  analyzed  Is  within  the  linear  lift  range  and  that  the 
panels  are  In  the  planar  (or  +)  orientation.  It  has  been  customarily  used  at 
angles  of  attack  where  nonlinear  lift  effects  occur.  But  with  a  significant 
overprediction  of  the  total  configuration  characteristics.  Hence,  It  is  not 
completely  suitable  for  missile  applications  where  the  angle  of  attack  regime 
required  Is  above  20°. 

An  alternative  approach,  and  that  selected  for  use  In  Missile  DATCOM.  Is 
the  "equivalent  angle  of  attack  method"  developed  by  Nielsen,  et  al  (Reference 
43).  This  method  assumes  that  the  panel  loading  for  a  given  panel  angle-of- 
attack  Is  unique.  With  this  method  the  panel  angle  of  attack  Is  computed  in¬ 
cluding  the  effect  of  panel  roll  orientation  with  respect  to  the  free  stream 
velocity  vector,  panel  proximity  to  the  fuselage  or  to  other  panels,  and  ex¬ 
ternal  vortex  flow  field  effects.  Then  the  Isolated  panel  characteristics  are 
interpolated  at  the  panel  quivalent  angle  of  attack  to  yield  the  panel  load 
when  mounted  on  a  body  In  combination  with  other  surfaces. 

This  concept  works  extremely  well.  It  Inherently  includes  the  nonlinear 
panel  loading  with  angle-of-attack,  from  the  isolated  panel  characteristics 
and  makes  the  synthesis  of  the  complete  configuration  aerodynamics  an  easy  to 
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apply  process.  The  difficult  task  Is  the  computation  of  the  panel  equivalent 
angle-of -attack,  but  this  Is  greatly  simplified  since  most  of  the  methodology 
Is  directly  derived  from  the  original  method  described  above. 


Figures  50  and  51  Illustrate  the  relative  Importance  of  component  Inter¬ 
ference  on  total  configuration  aerodynamics.  These  figures  demonstrate  that 
the  body,  lifting  surfaces,  and  carryover  Interference  effects  (neglecting  the 
Influence  of  vortices)  tend  to  overpredict  configuration  normal  force  by  15 
percent.  The  effect  of  body  and  wing  vortices  on  the  tail  surfaces  are  to 
reduce  total  normal  force,  since  they  decrease  the  tall  surface  local  angle- 
of-attack. 


Figure  52  summarizes  the  methods  selected  for  configuration  synthesis  of 
missile  configurations.  Note  that  only  the  body  vortex  strength  and  tracking 
methods  are  empirically  based;  this  Is  due  to  the  complexity  of  this  phenome¬ 
non  and  the  ready  availability  of  data  correlations  for  a  wide  range  of  mis¬ 
sile  body  shapes.  Table  4  summarizes  the  method  differences  between  those 
selected  from  the  Feasibility  Study  and  those  eventually  Incorporated  Into  the 
methods  compendium,  and  Figure  53  Illustrates  the  applicability  of  the  compo¬ 
nent  build  up  methods.  The  carryover  Interference  method  originally  selected 
were  those  documented  In  NACA  1307  (Reference  44);  they  were  replaced  by  the 
AIAA  Journal  methods  (References  45-47)  since  they  are  (1)  closed  form  solu¬ 
tions  as  a  function  of  afterbody  length  (which  was  not  Included  In  the  origi¬ 
nal  reference)  and  (2)  adjusted  to  account  for  panel  tralling-edge  sweep 
(which  the  original  method  did  not  address).  The  panel-panel  Interference 
correlation  developed  by  Nielsen  was  retained  for  use,  but  Is  limited  to 


planar  or  cruciform  fin  panels;  it  Is  not  applicable  to  multi-paneled  shapes 
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method  proposed  by  Darling  (Reference  48)  Is  not  approplate  for  the  "equiva¬ 
lent  angle  of  attack  method"  and  was  not  Incorporated.  Instead  Nielsen's 
method  has  been  extended  to  fin  sets  with  three  or  more  fins  by  mapping  the 
fins  adjacent  to  the  panel  being  analyzed  into  a  four  fin  geometry. 
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AMOtS  Of  ATTACK 


FIGURE  50  RELATIVE  IMPORTANCE  OF  THE  WING  VORTEX  CONTRIBUTION  TO  A 
CRUCIFORM  "+"  MISSILE  CONFIGURATION 
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FIGURE  51  RELATIVE  IMPORTANCE  OF  THE  WING  VORTEX  CONTRIBUTION  TO  A 


CRUCIFORM  "X"  MISSILE  CONFIGURATION 
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FIGURE  52  CONFIGURATION  SYNTHESIS  METHODS 
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TABLE  4  SYNTHESIS  METHODS 


Configuration  Synthwis 


1 1 • 4790 
7-20-84 


\  APPLICATION  SUITABLE 
?  APPLICATION  QUESTIONABLE, LIMITED 


FIGURE  53  COMPONENT  BUILD-UP  APPLICABILITY 


No  methods  were  selected  for  the  effect  of  fin  gaps  since  those  Influences 
have  been  demonstrated  to  be  small.  Some  available  experimental  data  indi¬ 
cates  that  the  resulting  error  is  less  than  10  percent  of  the  configuration 
normal  force  coefficient.  There  are  currently  no  methods  that  enable  the 
culatlon  of  this  effect  for  a  wide  range  of  fin  panels;  the  limited  empirical 
data  correlations  cannot  be  extended  to  handle  a  wide  variety  of  fin-body  com¬ 
binations  that  one  would  experience  in  the  missile  design  environment.  There¬ 
fore,  no  method  Is  included  to  predict  the  effect  of  fin  panel  "unporting" 
although  a  methods  development  task  should  undertaken  to  correct  this  de¬ 
ficiency. 


The  accuracy  of  the  methods  implement  id  were  demonstrated  by  comparison  to 
the  configurations  shown  in  Figure  54.  As  shown  in  Figires  55  and  56,  the 
predictive  accuracy  of  the  methods  is  excellent.  Note  that  predictive  accu¬ 
racy  was  determined  by  comparing  the  prediction  to  test  data  and  then  applying 
the  accuracy  criteria  defined  in  Figure  5. 
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FIGURE  54  AXISYMMETRIC  BODY  CONFIGURATIONS  INVESTIGATED 


6.1  EQUIVALENT  ANGLE  OF  ATTf.CK 

The  advantage  of  the  equivalent  angle-of-attuck  method  is  that  the  aero¬ 
dynamics  for  Isolated  fin  panels,  Including  any  nonllnearltles,  are  used  in 
the  prediction  of  total  configuration  aerodynamics.  The  Isolated  fin  panel 
data  Is  Interpolated  using  the  fin  local,  or  equivalent  angle-of-attack.  This 
angle-of -attack  Includes  the  effect  of  fin  d 'flection,  presence  of  a  body, 
presence  of  adjacent  fin  panels,  fin  roll  angU*  position  about  the  body,  and 
any  external  vortices.  The  accuracy  of  this  method  is  Illustrated  in  Figure 
57.  The  panel  loads  for  a  fin  are  plotted  versus  roll  angle  (measured  from 
top  vertical  enter)  at  two  angles  of  attack.  These  plots  show  excellent  pre¬ 
dictions  even  In  the  transonic  regime  using  the  equivalent  angle  of  attack 
method . 
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78 

91% 
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14 

71% 
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6 

83% 
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3 

100% 

AIM— 9L 

10 

100% 

LFRED 

43 

69% 

MAVERICK 

7 

43% 
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84 

95% 
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86% 
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356 

86% 

FIGURE  55  LONGITUDINAL  AERODYNAMICS  PREDICTIVE  ACCURACY 
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FIGURE  56  COMPARISON  OF  WING  IN  PRESENCE  OF  THE  BODY 
PREDICTIONS  WITH  EXPERIMENT  AT  MACH  3.36 
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FIGURE  57  TRANSONIC  PANEL  LOAD  PREDICTIONS  GOOD 


6.2  CARRYOVER  INTERFERENCE 


Carryover  interference  accounts  for  the  change  in  panel  angle  of  attack 
due  to  the  presence  of  a  body,  since  the  body  produces  a  flow  field  which 
varies  across  the  fin  span.  The  method  selected  for  Missile  DATCOM  are  those 
developed  In  NACA  1307  (Reference  44).  They  have  been  modified  to  reflect 
empirical  corrections  due  to  fin  traillng-edge  sweep,  since  the  NACA  1307 
methods  are  for  nonswept  traillng-edge  panels.  The  carryover  Interference 
due  to  angle-of-attack  method  employed  at  supersonic  speeds  are  those  develop¬ 
ed  by  Fan  and  Vlra,  Reference  46.  They  include  the  effect  of  fin  afterbody 
length.  The  methods  are  a  collection  of  closed-form  equations  which  reproduce 
the  NACA  1307  results  for  no  afterbody  and  full  afterbody. 

The  effect  of  panel  incidence  on  carryover  interference  are  exactly  those 
developed  in  NACA  1307.  No  other  method  was  found  that  enables  the  prediction 
of  this  effect  other  than  empirical  data  correlations. 

The  effect  of  roll  angle  (combined  pitch  and  yaw  angle-of-attack)  are  com¬ 
puted  using  the  method  developed  by  Nielsen  (Reference  43)  from  slender  body 
theory.  Unfortunately,  the  method  addresses  only  planar  or  cruciform  fin 
panels.  Although  one  could  Interpolate  between  the  planar  and  cruciform  re¬ 
sults  to  model  equally  spaced  triform  fin  arrangements,  other  fin  combinations 
are  not  addressed.  A  method  should  be  developed  to  accommodate  these  types  of 
general  fin  arrangements. 

There  Is  need  for  further  method  development  of  an  angle-of-attack  correc¬ 
tion  to  body-wing  carry-over  for  geometries  where  the  body  Is  dominant  (Refer¬ 
ence  50).  This  is  Illustrated  In  Figure  57  which  shows  normal  force  coeffi¬ 
cient  In  presence  of  the  body  versus  angle-of-attack  for  various  aspect 
ratios  and  body  radius  to  span  ratios  (Reference  50).  This  figure  illustrates 
that  the  error  increases  with  a  decrease  in  aspect  ratio,  and  decreases  with  a 
decrease  in  body  radius  to  span  ratio. 
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6.3  VORTEX  TRACKING  AND  STRENGTH 


Two  vortex  sources  are  Included  In  the  Missile  Datcom  analysis  of  the  com¬ 
plete  configuration  aerodynamics.  These  are  (1)  vortices  which  are  shed  from 
the  body  nose,  and  (2)  vortices  that  are  shed  from  any  forward  lifting  sur¬ 
faces.  Methods  for  vortices  which  originate  from  the  afterbody  have  not  been 
Included.  This  was  done  to  simplify  the  analysis  since  complex  flow  calcula¬ 
tions  are  required  to  adequately  model  these  vortices.  The  Influence  of  these 
vortices  on  aft  mounted  fin  panels  are  minor  and  can  be  Ignored  In  preliminary 
missile  design.  It  is  hoped  that  an  easy  to  apply  technique  will  be  develop¬ 
ed  to  completely  model  the  vortices  remaining  In  the  flow  field. 

6.3.1  Body  Vortices 

The  body  vortices  are  tracked  using  a  method  developed  by  Nielsen,  et.al. 
(Reference  51).  It  relies  upon  an  extensive  empirical  correlation,  and  has 
been  observed  to  be  adequate  for  preliminary  design.  The  strength  of  the 
vortices  are  computed  using  the  lift  generated  by  the  body. 

6.3.2  Lifting  Surface  Vortices 

The  fin  vortices  are  computed  using  the  methods  outlined  in  the  U.S.  Air 
Force  DATCOM.  The  vortices  are  tracked  assumming  that  they  follow  the  path  of 
the  free  stream  velocity  vector.  This  technique  showed  good  predictions  ex¬ 
cept  when  the  vortex  core  generated  by  a  forward  fin  impinged  upon  an  aft 
fin.  A  new  model  was  developed  which  established  a  vortex  core  radius  of  .255 
of  the  fin  span  generating  the  vortex.  The  vortex  is  then  tracked  along  the 
free  stream  unless  it  Impinges  on  a  trailing  fin  In  which  case  it  Is  offset  by 
the  core  radius.  Figure  58  shows  the  improvement  of  this  method  over  the 
original  method.  The  vortex  strength  is  computed  from  the  fin  lift.  Note 
that  fin  lift  Is  dependent  upon  panel  local  angle  of  attack,  which  In  turn  Is 
dependent  upon  any  external  vortices.  Hence,  In  order  to  accurately  compute 
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pane]  lift,  and,  vortex  strength,  the  pane]  ]oads  for  each  panel  on  the  con¬ 
figuration  must  be  calculated  starting  from  the  nose  of  the  vehicle  and  work¬ 
ing  aft.  Therefore,  the  effect  of  vortices  originating  ahead  of  a  fin  set  Is 
Included  In  the  calculation  of  that  fin  set's  equivalent  angle  of  attack. 


GOfllLLION  B-W-T,  'X'  CONFIGURATION 


FIGURE  58  EFFECT  OF  NEW  VORtEX  CORE  MODEL 


7.0  DYNAMIC  DERIVATIVE  METHOD  SELECTION 


This  section  summarizes  the  methods  selected  for  dynamic  derivatives.  It 
Is  divided  into  three  parts:  body  alone  methods,  fin  alone  methods  and  con¬ 
figuration  synthesis.  Each  section  is  further  subdivided  by  coefficient. 

The  available  methodologies  are  for  conventional  missile  shapes  at  low 
angle  of  attack.  The  methods  were  not  as  numerous  as  those  for  the  static 
coefficients,  and  no  methods  were  found  that  specifically  address  the  effect 
of  body  Incidence.  Longitudinal  dynamic  derivative  methods  for  both  body-  and 
fin-alone,  and  Magnus  coefficient  methods  for  body  alone  were  selected. 
Figure  59  summarizes  the  methods  selected. 

7.1  BODY  ALONE 

The  available  body-alone  methodology  suitable  for  preliminary  design  is 
derived  from  slender  body  theory,  or  correlations  of  large  amounts  of  test 
data.  The  large  effects  of  viscosity  over  slender  bodies  at  anything  but 
small  angles  of  attack  make  slender  body  results  applicable  only  in  the  linear 
angle-of-attack  range.  Although  viscous  effects  such  as  boundary  layer  tran¬ 
sition  have  been  shown  to  have  significant  effects  on  dynamic  stability,  the 
mathematical  complexity  involved  in  these  solutions  renders  them  unsuitable 
for  preliminary  design. 

7.1.1  Normal  Force  Due  to  Pitch  Rate 

The  method  selected  Is  that  of  the  United  States  Air  Force  Stability  and 
Control  uATCOM.  The  method  used  subsonically,  transonical ly  and  super¬ 
sonically  were  derived  from  slender  body  theory.  Also  there  is  a  hypersonic 
method  derived  from  simple  Newtonian  theory. 
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SUMMARY  OF  DYNAMIC  DERIVATIVE  METHODOLOGY 


7.1.2  Pitching  Moment  Due  to  Pitch  Rate 


As  recommended  in  the  Missile  DATCOM  Feasibility  Report,  the  method  of 

Ericsson  was  chosen  (Reference  54).  Ericsson  extended  the  simple  slender  body 

theory  and  Newtonian  theory  to  more  properly  account  for  Mach  number.  The 

method  computes  the  total  damping  in  pitch  derivative  C  +  C  . .  These 

m,  ma 

two  are  computed  as  their  sum  In  both  the  handbook  and  the  computer  code. 
Figure  60  shows  good  comparisons  for  body  alone  pitch  damping.  For  the  hand¬ 
book  a  less  mathematically  Intensive  method  than  that  of  Reference  55  was 

chosen  to  compute  C  supersonically.  The  method  is  that  of  DATCOM  for 

nu 

supersonic  body  alone  C  .  If  C  alone  and  not  the  total  pitch  damping 

m*  mo 

derivative  is  desired  the ’  method  of  DATCOM  may  be  used  in  all  Mach  regimes. 
For  hypersonic  Mach  numbers,  the  method  of  DATCOM  was  selected  which  applies 
simple  Newtonian  theory  to  calculate  C  of  cones  with  or  without  spherical 

nIa 

nose  caps  and  conic  frustums.  ' 

7.1.3  Normal  Force  Due  to  Rate  of  Change  of  Angle  of  Attack 

The  method  chosen  Is  that  of  DATCOM  which  was  derived  using  slender  body 
theory.  This  method  is  applicable  to  all  Mach  regimes. 

7.1.4  Pitching  Moment  Due  to  Rate  of  Change  of  Angle  of  Attack 

As  mentioned  In  Section  7,1.2,  the  method  selected  Is  that  of  Ericsson 
which  gives  the  total  damping  In  pitch  derivative  C  +  C  ..  For  the 

m  a  life 

handbook,  a  different  method  was  selected  for  Mach  numbers  gr  iter  than  one, 
which  Is  better  suited  for  hand  calculations.  The  method  is  that  of  DATCOM 
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desired  and  not  the  total  pitch 
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damping  derivative  the  method  of  DATCOM  may  be  used  in  all  Mach  regimes. 


7.1.5  Side  force  Due  to  Roll  Rate 


An  empirical  method  was  selected  from  Reference  56.  The  method  is  derived 
from  wind  tunnel  and  ballistic  range  test  on  spin  stabilized  projectiles. 
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FIGURE  60  BODY  ALONE  PITCH  DAMPING  COMPARISON 


This  method  Is  good  for  bodies  whose  length  is  between  2.5  and  10  calibers, 
whose  nose  length  Is  between  1.2  to  5.5  calibers  and  whose  boattall  length  Is 
between  0  to  1  caliber. 

7.1.6  Yawing  Moment  Due  to  Roll  Rate 

The  method  selected,  which  gives  C  /sina,  Is  an  empirical  method 
based  on  wind  tunnel  and  ballistic  range  test  of  spin  stabilized  projectiles 
(Reference  56).  Caution  should  be  taken  using  this  method  since  C^.sina 
versus  slna  does  not  continue  to  behave  In  polynomial  fashion  at  higher 
angles  of  attack  (see  Figure  61).  This  method  should  not  be  applied  beyond 
the  second  positive  maximum  (counting  the  one  at  a-0  as  the  first).  Addi¬ 
tional  limitations  on  this  method  are  listed  In  Section  7.1.5. 

7.1.7  Rolling  Moment  Due  to  Roll  Rate 

The  roll  damping  derivative  method  selected  Is  that  of  Reference  57.  The 
method  Is  empirical  and  was  derived  from  correlations  of  wind  tunnel  and 
ballistic  range  test.  Limitations  on  this  method  are  listed  In  Section  7.1.5. 


FIGURE  61  LIMIT  OF  YAWING  MOMENT  DUE  TO  ROLL  RATE  METHOD 


7.2  FIN  ALONE 


The  only  available  methods  suitable  for  use  In  Missile  OATCOM  to  compute 
fin  alone  dynamic  derivatives  are  the  methods  of  the  United  States  Air  Force 
Stability  and  Control  OATCOM.  The  methods  are  based  on  lifting  surface  theory 
for  subsonic  speeds  and  on  linearized  theory  for  supersonic  speeds.  The 
methods  are  thus  limited  to  conditions  for  which  the  flow  Is  attached  over  the 
surface  of  the  fin,  l.e.,  the  linear  lift  range.  These  methods  were  developed 
for  airplanes  which  generally  have  larger  aspect  ratio  lifting  surfaces. 
There  are  some  limits  on  AR  (where  Is  the  Mach  similarity  factor)  for  the 
supersonic  methods.  In  the  supersonic  regime  the  methods  were  developed  for 
beta  times  aspect  ratio  greater  than  four.  For  typical  missile  aspect  ratios 
at  low  supersonic  Mach  numbers  beta  times  aspect  ratio  Is  commonly  less  than 
four.  For  sample  cases,  however,  the  trends  beyond  this  limit  look  reason¬ 
able.  However  the  accuracy  of  this  method  for  lower  aspect  ratios  cannot  be 
determined  because  only  a  small  amount  of  experimental  data  Is  available,  and 
that  data  Is  for  subsonic  Mach  numbers. 

7.2.1  Normal  Force  and  Pitching  Moment  Due  to  Pitch  Rate 

The  methods  selected  are  from  OATCOM.  The  method  used  depends  upon  Mach 
regime  and  planform  geometry.  In  the  transonic  regime  a  linear  fairing  Is 
used  between  the  value  of  the  coefficient  at  Mach  equal  to  0.8  and  1.2. 
Supersonically,  when  the  leading  edge  Is  subsonic,  there  are  methods  for  fins 
with  taper  ratio  of  0  or  taper  ration  between  0.25  and  1,  there  are  no  methods 
for  taper  ratio  between  0  and  0.25.  In  the  computer  code  a  linear  fairing  Is 
used  between  the  value  of  the  coefficient  at  a  taper  ratio  of  0  and  of  0.25. 
There  are  separate  methods  for  subsonic  and  supersonic  leading  edges.  During 
the  method  change  over  when  the  leading  edge  is  sonic  or  slightly  higher  the 
methods  do  not  match,  and  the  supersonic  leading  edge  method  yields  unrealis¬ 
tic  results. 
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7.2.2  Normal  Fore*  and  Pitching  Moment  Due  to  Rate  of  Change  of  Angle  of 
Attack 

The  methods  selected  are  from  DATCOM.  The  subsonic  method  is  applicable 
to  fins  with  triangular  planforms.  Supersonically,  when  the  fin  has  a  sub¬ 
sonic  leading  edge,  there  are  no  methods  for  taper  ratio  between  0  and  0.25; 
methods  are  available  for  a  taper  ratio  of  0  and  taper  ratios  between  0.25  and 
1.  Supersonically  when  the  fin  has  l  supersonic  leading  edge,  the  computer 
code  uses  a  table  look-up  on  design  charts  to  compute  these  coefficients. 

7.3  CONFIGURATION  SYNTHESIS 

Configuration  synthesis  is  done  by  determining  the  contribution  of  each 
component  using  the  body  alone  and  fin-alone  methods  previously  discussed. 
The  contribution  of  each  fin  to  a  fin  set  is  determined  by  taking  Into  account 
the  fin's  position  on  the  body.  The  contribution  of  each  fin  Is  then  summed 
to  compute  the  contribution  of  a  fin  set.  Body  on  fin  and  fin  on  body  Inter¬ 
ference  factors  are  computed  using  slender  body  theory.  Then  the  coefficient 
for  the  total  configuration  Is  the  body  contribution  plus  the  product  of  the 
Interference  factor  and  the  fin  set  contribution. 

Figures  62  and  63  show  comparisons  of  two  configurations.  The  Sidewinder 
comparisons  (Figure  62)  for  the  body-canard-tall  In  *+■  configuration  show 
errors  between  19  and  27  percent.  Both  the  body  tall  configurations  (Flnner 
and  Sidewinder)  over  predict  damping  n  pitch.  The  error  for  the  sidewinder 
body-tall  configuration  ranges  from  8.6  to  20  percent. 
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E  63  CONE-CYLINDER  TAIL  PITCH  DAMPING  COMPARISON 


8.0  COMPUTER  PROGRAM 


The  Missile  DATCOM  computer  code  was  written  In  ANSI  FORTRAN  using  struc¬ 
tured  programming  techniques  Including  In-code  documentation.  Input  to  the 
code  Is  accomplished  by  namelist  Inputs  and  control  cards  that  are  read  by  a 
namelist  emulator  which  Is  part  of  the  code.  The  namelist,  emulator  was 
developed  because  not  all  machines  have  FORTRAN  namelist  and  the  code  Is  not 
easily  converted  to  fixed  field  Inputs.  Once  read,  the  inputs  are  checked  for 
valid  variables  and  major  errors.  If  an  Input  error  Is  found,  a  diagnostic 
message  Is  written  to  the  output  file  describing  what  the  error  Is  and  where 
it  Is. 

The  computer  code  can  be  run  either  in  "un-overlayed"  or  In  "overlayed* 
mode.  In  " un-over layed"  mode,  the  required  core  memory  is  approximately 
277,000  octal  words  on  CDC  machines.  If  the  program  Is  "overlayed"  the  core 
requirement  on  CDC  machines  Is  76,000  octal  words.  On  the  virtual  memory  VAX 
machine  the  core  requirement  Is  735,744  bytes.  The  code  can  be  easily  con¬ 
verted  from  FORTRAN  IV  to  FORTRAN  V.  The  user  is  referred  to  Volume  2  of  the 
Program  User's  Guide  for  more  details. 

The  computer  code  was  developed  In  parallel  with  the  handbooks.  Most  of 
the  methods  in  the  handbooks  are  incorporated  Into  the  computer  code.  In  some 
Instances  a  more  accurate  but  mathematically  more -Intensive  method  was  In¬ 
corporated  Into  the  code. 

The  methods  selected  were  coded  Into  Individual  subroutines  cr  In  the  case 
of  complex  methods  Into  a  group  of  subroutines  called  method  modules.  To  re¬ 
place  a  method,  a  new  method  Is  coded  Into  a  subroutine  and  inserted  in  place 
of  the  old  subroutine.  In  the  case  of  method  modules  the  substltitlon  Is  a 
little  more  complex  but  is  still  easily  accomplished.  Because  the  code  was 
developed  using  top-down  design,  most  of  the  control  logic  is  In  the  upper 
level  of  the  code  and  would  not  require  a  change. 
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The  computer  code  also  has  two  additional  capabilities  which  allow  for 
more  detailed  analysis  of  a  configuration,  experimental  data  substitution  and 
configuration  Incrementing. 

Experimental  data  substitution  allows  the  user  to  Input  experimental  data 
for  any  part  of  the  configuration  (body,  fin  set  1,  etc.)  or  any  partial  con¬ 
figuration  (body  +  1  fin  set,  etc.).  For  example  the  user  may  Input  body 
alone  experimental  data  and/or  body  +  1  fin  set  experimental  data  for  a  con¬ 
figuration  having  a  body  +  2  fin  sets.  The  experimental  data  is  substituted 
into  the  appropriate  common  block  and  replaces  the  computed  theoretical  co¬ 
efficients.  The  substutlted  coefficients  are  then  used  In  configuration  syn¬ 
thesis. 

Configuration  incrementing  allows  the  user  to  input  experimental  data  for 
a  configuration  In  the  first  case  of  a  run  and  then  vary  so<ae  components  of 
the  configuration.  In  subsequent  cases,  correction  factors  are  applied  to  the 
static  aerodynamic  coefficients.  When  incrementing  all  cases  of  a  run  must  be 
for  the  same  configuration  (l.e.,  If  the  first  case  Is  a  body  +  1  fin  set  then 
the  following  case  must  also  be  a  body  +  1  tin  set).  In  the  first  case,  ex¬ 
perimental  data  Is  substituted  for  the  configuration  being  run.  Comparisons 
are  then  made  between  the  experimental  and  theoretical  coefficients,  and 
correction  terms  are  calculated.  In  subsequent  cases  of  that  run,  the  correc¬ 
tion  factors  are  used  to  Increment  the  theoretical  coefficients  of  the  con¬ 
figuration. 
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